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Abstract: We use the antenna subtraction method to isolate the mixed real- virtual 
infrared singularities present in gluonic scattering amplitudes at next-to-next-to-leading 
order. In a previous paper, we derived the subtraction term that rendered the double 
real radiation tree-level process finite in the single and double unresolved regions of phase 
space. Here, we show how to construct the real-virtual subtraction term using antenna 
functions with both initial- and final-state partons which removes the explicit infrared 
poles present in the one-loop amplitude, as well as the implicit singularities that occur in 
the soft and collinear limits. As an explicit example, we write down the subtraction term 
that describes the single unresolved contributions from the five-gluon one-loop process. 
The infrared poles are explicitly and locally cancelled in all regions of phase space prior 
to integration, leaving a finite remainder that can be safely evaluated numerically in 
four-dimensions. We show numerically that the subtraction term correctly approximates 
the matrix elements in the various single unresolved configurations. 
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1. Introduction 



In hadronic collisions, the most basic form of the strong interaction at short distances 
is the scattering of a coloured parton off another coloured parton. Experimentally, such 
scattering can be observed via the production of one or more jets of hadrons with large 
transverse energy. In QCD, the (renormalised and mass factorised) inclusive cross section 
has the form, 

where the probability of finding a parton of type i in the proton carrying a momentum 
fraction £ is described by the parton distribution function fi{£,,np)d£ and the parton- level 
scattering cross section d&ij for parton i to scatter off parton j normalised to the hadron- 
hadron flux 1 is summed over the possible parton types i and j. As usual and fip are 
the renormalisation and factorisation scales. 

The infrared-finite partonic cross section for a parton of type i scattering off parton of 
type h has the perturbative expansion 

d*« = daff + (^M) da^° + (^=M) 2 da™° + Otf) (1.1) 

where the next-to-leading order (NLO) and next-to-next-to-leading order (NNLO) strong 
corrections are identified. Note that the leading order cross section is 0{o? s ). 

The single-jet inclusive and dijet cross sections have been studied at NLO [1-6] and 
successfully compared with data from the high energy frontier at the TEVATRON [7-9] 
and at the LHC [10, 11]. A particular success is the determination of the strong coupling 
constant as a function of the jet transverse energy [12-14]. 

The theoretical prediction may be improved by including the NNLO perturbative pre- 
dictions. This has the effect of (a) reducing the renormalisation scale dependence and (b) 
improving the matching of the parton level theoretical jet algorithm with the hadron level 
experimental jet algorithm because the jet structure can be modelled by the presence of a 
third parton. The resulting theoretical uncertainty at NNLO is estimated to be at the few 
per-cent level [15]. 

Any calculation of these higher-order corrections requires a systematic procedure for 
extracting the infrared singularities that arise when one or more final state particles become 
soft and/or collinear. These singularities are present in the real radiation contribution at 
NLO, and in double real radiation and mixed real- virtual contributions at NNLO. 

Subtraction schemes are a well established solution to this problem and several meth- 
ods for systematically constructing general subtraction terms have been proposed in the 
literature at NLO [16-20] and at NNLO [21-37]. Another NNLO subtraction type scheme 
has been proposed in [38]. It is not a general subtraction scheme, but can nevertheless deal 

lr rhe partonic cross section normalised to the parton-parton flux is obtained by absorbing the inverse 
factors of £1 and £2 into d&ij. 
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with an entire class of processes, those without coloured final states, in hadron-hadron colli- 
sions. It has been explicitly applied to several observables [38-43]. In addition, there is the 
completely independent sector decomposition approach which avoids the need for analyti- 
cal integration and which has been developed for virtual [44-47] and real radiation [47-51] 
corrections to NNLO, and applied to several observables already [52-55]. 

We will follow the NNLO antenna subtraction method which has been fully developed 
for the production of massless partons in electron-positron annihilation [21,56,57], leading 
to the successful description of the infrared structure of three-jet events at NNLO [58, 
59] and the subsequent numerical calculation of the NNLO corrections to event shape 
distributions [60-63], the moments of event shapes [64,65] and jet rates [66,67]. 

The formalism has been extended to include massive fermions at NLO [68, 69] and 
recently Bernreuther et al applied the antenna subtraction scheme to electron-positron 
annihilation with massive final state quarks at NNLO [70]. 

Here, we are more interested in the application of the antenna subtraction scheme 
to hadron-hadron collisions, and specifically to processes such as pp — > jet + X, pp — > 
V + jet + X (with V = or Z) and pp — > H + jet + X. The formalism has been fully 
worked out for initial state hadrons at NLO [71] and further developed at NNLO [72-74]. 
Within this method, the subtraction terms are constructed from so-called antenna functions 
which describe all unresolved partonic radiation (soft and collinear) between a hard pair 
of radiator partons. The hard radiators may be in the initial or in the final state, and 
in the most general case, final-final (FF), initial-final (IF) and initial- initial (II) antennae 
need to be considered. Within the antenna subtraction method, the subtraction terms and 
therefore the antennae also need to be integrated over the unresolved phase space, which is 
different in the three configurations. All of the integrals for the final-final and initial-final 
configurations are fully known. Work is presently underway to determine the few remaining 
integrals of tree-level four parton antenna for the initial-initial configuration. A first step 
in that direction has been realised in [73]. Once completed, the computation of the NNLO 
corrections for jet production at the LHC is in range. 

Suppressing the labels of the partons in the initial state of the hard scattering, the 
general form for the subtraction terms for an m-particle final state at NNLO is given 
by [21]: 



da NNLO = / (da§NLO ~ d&NNLo) + / &°NNLO 

+ / fa NNLO ~ d °NNLo) + / d °NNLO + / d<5 " 



+ / d °NNLO + j &°NNLO- ( L2 ) 

Here, da NNLO denotes the subtraction term for the (m+2)-parton final state which behaves 
like the double real radiation contribution dcrjy^y LO in all singular limits. Likewise, d &NNLO 
is the one-loop virtual subtraction term coinciding with the one-loop (m + l)-final state 
dcrjvNLO in all singular limits. The two-loop correction to the m-parton final state is 
denoted by da]^ L0 . In addition, when there are partons in the initial state, there are two 
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mass factorisation contributions, da NNLO and da NNLO , for the (m + 1)- and m-particle 
final states respectively. 

In order to construct a numerical implementation of the NNLO cross section, the 
various contributions must be reorganised according to the number of final state particles, 

da nnlo = / [da nnlo ~ da NNLO ] 



+ / [d&NNLO ~ daJfNLo] 

Jd<S> m+1 

+ / [d^NNLO ~ da NNLO ] , (1.3) 

where the terms in each of the square brackets is finite and well behaved in the infrared 
singular regions. More precisely, 

d&NNLO = da N S NLO - J^d&^NLO ~ d a NNLO> 0-^) 

da nnlo = ~~ J^daYfNLO ~ J da^ NLO - da^^'l Q . (1.5) 

Note that because integration of the subtraction term da NNLO gives contributions to both 
the (m+1)- and m-parton final states, we have explicitly decomposed the integrated double 
real subtraction term into a piece that is integrated over one unresolved particle phase space 
and a piece that is integrated over the phase space of two unresolved particles, 

/ da NNLO = / da^ NLO + / / da S / NLQ . (1.6) 

In a previous paper [75], the subtraction term da NNLO corresponding to the leading 
colour pure gluon contribution to dijet production at hadron colliders was derived. The 
subtraction term was shown to reproduce the singular behaviour present in d<r|^y i0 in all 
of the single and double unresolved limits. 

It is the purpose of this paper to construct the appropriate subtraction term da NNLO to 
render the leading colour five-gluon contribution da 1 ^^ explicitly finite and numerically 
well behaved in all single unresolved limits. Our paper is organised in the following way. In 
Section 2, the general structure of f ± da^NLO anc ^ ^nnlo are discussed and analysed. The 
coupling constant renormalisation to remove the UV singularities and the mass factorisation 
of the initial-state singularities into the parton distributions are described in Section 3. In 
Section 4 we turn our attention to the specific process of gluon scattering at NNLO. Our 
notation for gluonic amplitudes is summarised in Section 4.1 and the one- loop five-gluon 
amplitudes discussed in Section 4.2. 

There are two separate configurations relevant for gg — > ggg scattering depending on 
whether the two initial state gluons are colour-connected or not. We denote the configura- 
tion where the two initial state gluons are colour-connected (i.e. adjacent) by IIFFF, while 
the configuration where the colour ordering allows one final state gluon to be sandwiched 
between the initial state gluons is denoted by IFIFF. Explicit forms for the real-virtual 
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subtraction term dajj NLO for these two configurations are given in Section 5 where the 
cancellation of explicit poles is made manifest. The validity of the subtraction term is 
tested numerically in Section 6 by studying the subtracted one-loop matrix elements in all 
of the single unresolved limits. In particular, we show that in all cases the ratio of the 
finite part of the real- virtual and subtraction terms approaches unity. Finally, our findings 
are summarised in Section 7. 

Four appendices are also enclosed. Appendix A summarises the phase space map- 
pings for the final-final, initial-final and initial-initial configurations. Appendix B gives 
a description of the tree-level antennae appearing in the subtraction terms in both their 
unintegrated form (Appendix B.l) and after integration over the unresolved phase space 
(Appendix B.2). The unintegrated one-loop antenna is given in Appendix B.3. Appendix C 
contains a modified form for the wide angle soft subtraction terms present in daf iNLO , while 
formulae relating to the mass factorisation contribution are given in Appendix D. 



2. Real- Virtual antenna subtraction at NNLO 



In this paper, we focus on the kinematical situation of the scattering of two massless 
coloured partons to produce massless coloured partons, and particularly the production 
of jets from gluon scattering in hadronic collisions. To establish some notation, consider 
the leading-order parton-level contribution from the (m + 2)-parton processes to the m-jet 
cross section at LO in pp collisions, 

pp — >■ m jets (2.1) 

which is given by 

d& LO = MlO ^ d ®rn(P3, ■ ■ ■ , Pm+2] Pi , P2) J~ 

O 777 

perms 

x \M m +2(l 2 . . . , m + 2)\ 2 jW(p 3 , . . . , Pm+2 ). (2.2) 

To make the subsequent discussion more general, we denote a generic tree-level (m + 2)- 
parton colour ordered amplitude by the symbol A1 m +2(1, 2 . . . , m + 2), where 1 and 2 
denote the initial state partons of momenta p\ and p 2 while the m-momenta in the final 
state are labeled p$, ■ ■ ■ ,p m +2- For convenience, and where the order of momenta does not 
matter, we will often denote the set of (m + 2)-momenta {pi, . . . ,p m +2} by {p} m+ 2- The 
symmetry factor S m accounts for the production of identical particles in the final-state. 

At leading colour, the colour summed squared matrix elements are determined by the 
squares of the individual colour ordered amplitudes where the sum runs over the various 
colour ordered amplitudes. 2 For gluonic amplitudes this is the sum over the group of 
non-cyclic permutations of n symbols. The normalisation factor Nlo includes the hadron- 
hadron flux-factor, spin and colour summed and averaging factors as well as the dependence 
on the renormalised QCD coupling constant a s . 

2 The sub-leading colour contributions are implicitly contained in the squared matrix elements. 
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d<J> m is the 2 — > m particle phase space: 

d$m{P3, ■ ■ ■ ,Pm+2;Pl,P2) = 

d d_1 P3 d d-1 p m+2 . d d , . . . 

2^(2^ • • • 2W2vr)^ (27r) ^ + " ** " ' ' ' " Pm+2) ' ^ 

The jet function Jm^({p} n +2) defines the procedure for building m-jets from n final state 
partons. Any initial state momenta in the set {p} n+ 2 are simply ignored by the jet algo- 

(n) 

rithm. The key property of Jm is that the jet observable is collinear and infrared safe. 

In a previous paper [75], we have discussed the NNLO contribution coming from pro- 
cesses where two additional partons are radiated - the double real contribution da^^ LO 
and its subtraction term d&ff NLQ . da N ^f LO involves the (m + 4)-parton process at tree 
level and is given by, 

dcr§N LO = J^NNLO Y] d$ m+2 (p 3 , . . . , p m+4 ; Pi, P2) 
perms 

x \M m+4 (l 2 . . . , m + 4)| 2 jt +2) (P3, Mm+4). (2-4) 



In this paper, we are concerned with the NNLO contribution coming from single radiation 
at one-loop, i.e. the (m + 3)-parton process, da^f L0 . 

In our notation, the one-loop (m+3)-parton contribution to m-jet final states at NNLO 
in hadron-hadron collisions is given by 

d&NNLO =NnNLO d$ m+l(P3, • • • , Pm+3] Pi, P2) 



perms 

x IMl^l, . . . , m + 3)| 2 Jt +1) ({ P } m+3 ) , (2.5) 



Sm+l 

,1 ™ , Q\|2 T (m+l)ff m \ _\ 

where we introduced a shorthand notation for the interference of one-loop and tree-amplitudes, 

|^ +3 (i,...,m + 3)| 2 = 2Re (^ +3 (i,...,m + 3)^* +3 (i,...,m + 3)) , (2.6) 

which explicitly captures the colour-ordering of the leading colour contributions. The sub- 
leading contributions in colour are implicitly included in (2.6) but will not be considered 
in detail in this paper. As usual, the sum is the appropriate combination of colour or- 
dered amplitudes. For gluonic amplitudes, this is the sum over the group of non-cyclic 
permutations of n symbols. 

The normalisation factor Nlo depends on the specific process and parton channel 
under consideration. Nevertheless, N^ NLO and ^nnlo are related both to each other and 
to Nlo f° r an Y number of jets and for any partonic process by 

KfRV _ K r (a s N\ 2 C(er 

Mnnlo - Nlo \~w ) W {) 

\fRR — \r (*sN\ 2 C(ef 

Mnnlo - ^lo I I c^\2> ^ ' 

Mnnlo = Mnnlo *C(e), (2-9) 
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where 

C(e) = (4vr) e ^, (2.10) 
C(e) = (4^) e e" £7 . (2.11) 

As expected, each power of the (bare) coupling is accompanied by a factor of C(e). In this 
paper, we are mainly concerned with the NNLO corrections to (2.1) when m = 2 and for 
the pure gluon channel. The normalisation factor Mlo will be given for this special case 
in Section 4. 

Eq. (2.5) contains two types of infrared singularities. The renormalised one-loop virtual 
correction A^^ +3 to the (m + 3)-parton matrix contains explicit infrared poles, which can 
be expressed using the infrared singularity operators defined in [76, 77]. On the other 
hand, the requirement of building m-jets from (m + l)-partons allows one of the final state 
partons to become unresolved, leading to implicit local infrared singularities which become 
explicit only after integration over the unresolved patch of the final state (m + l)-parton 
phase space. The single unresolved infrared singularity structure of one- and two-loop 
amplitudes has been studied in [78-90]. 

As discussed in Section 1, in order to carry out the numerical integration over the 
(m + l)-parton phase, weighted by the appropriate jet function, we have to construct an 
infrared subtraction term 3 dcr^ NLO which 

(a) removes the explicit infrared poles of the virtual one-loop (m + 3)-parton matrix 
element. 

(b) correctly describes the single unresolved limits of the virtual one-loop (m + 3)-parton 
matrix element. 

The subtraction term has three components; 

where J l do S ^ l NLO is derived from the double real radiation subtraction term dofj NLO in- 
tegrated over the phase space of one unresolved particle. Part of this contribution cancels 
the explicit poles in the virtual matrix element, while the real- virtual subtraction term 
d<7 jvwLO accounts for the single unresolved limits of the virtual matrix element. 

In the following subsections we shall present the general structure of J l da^ l NLO and 
daJfffLQ. The remaining poles are associated with the initial state collinear singularities 
and are absorbed by the mass factorisation counterterm da NN ^ which will be presented 
explicitly in Section 3. 

A key element of the antenna subtraction scheme is the factorisation of the matrix 
elements and phase space in the singular limits where one or more particles are unresolved. 



3 Strictly speaking, daJf NLO is not a subtraction term since it adds back part of the the double radiation 
subtraction term d&ff NLO integrated over the phase space of a single unresolved particle. Nevertheless, 
since it contains all the terms needed to render the (m + l)-particle final state finite, it is convenient to call 
it the real-virtual subtraction term. 
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In determining the various contributions to daJf NLO , we shall therefore specify the uninte- 
grated and/or integrated antennae and the reduced colour ordered matrix-element squared 
involved. The factorisation is guaranteed by the momentum mapping described in the Ap- 
pendix A. For conciseness, only the redefined hard radiator momenta will be specified in 
the functional dependence of the matrix element squared. The other momenta will simply 
be denoted by ellipsis. 

In order to combine the subtraction terms and real- virtual matrix elements, it is con- 
venient to slightly modify the phase space, such that 

dxi dx 2 



d °NNLO = NnNLO ^2 J d ®rn+l(P3,---,Pm+3]XlPl,X2P2)^r- 



perms 



S m+ i Xi X 2 



x IMl^i, . . . , m + 3)| 2 5(1 - Xl )S(l - x 2 ) J^ +1 \{p} m+3 ) . (2.13) 

The integration over x\ and x 2 reflects the fact that the subtraction terms contain contri- 
butions due to radiation from the initial state such that the parton momenta involved in 
the hard scattering carry only a fraction Xi of the incoming momenta. In general, there 
are three regions: the soft {x\ = x 2 = 1), collinear (x\ = 1, x 2 ^ 1 and x\ ^ 1, x 2 = 1) 
and hard (x\ ^ 1, x 2 ^ 1). The real- virtual matrix elements only contribute in the soft 
region, as indicated by the two delta functions. 

In Sections 2.1 and 2.2, we discuss the first two terms that contribute to d&Jf NLO 
given in Eq. (2.12), namely Jida^ l NLO and da^ff L0 . The final contribution da^^lo is 
discussed in Section 3. 

2.1 Contribution from integration of double real subtraction term: dafl]^ LO 

There are five different types of contributions to dafj NLO according to the colour connection 
of the unresolved partons: 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 

(c) Two unresolved partons that are not colour-connected but share a common radiator 
(almost colour-connected) . 

(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour- unconnected) . 

(e) Compensation terms for the over subtraction of large angle soft emission. 

Each type of contribution takes the form of antenna functions multiplied by colour ordered 
matrix elements. The various types of contributions are summarised in Table 1. We see 
that the a, c and e types of subtraction term, as well as the 6-type terms that are products 
of three-particle antenna, can be integrated over a single unresolved particle phase space 
and therefore contribute to the (m + l)-particle final state so that, 

j i d °NNLO = ^ d& NNLO + ^ d& NNLO + ^ d °NNLO- ( 2 - 14 ) 
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a 


b 


b, c 


d 


e 


ua NNLO 

II d °NNLO 

f A ~ S ' 2 
J2 aa NNLO 


X!\M° m+ 3\ 2 


x!\M° m+2 \ 2 

X2\M° n+2 \ 2 


X$X$\M° m+2 \-> 
^X°\Ml +2 \* 


XlX$\Ml +2 f 
X§X°\M° m+2 \ 2 


SX$\M a m+2 \ 2 
SX°\Ml +2 f 



Table 1: Type of contribution to the double real subtraction term daff NLO , together with the 
integrated form of each term. The unintegrated antenna and soft functions are denoted as X®, 
X® and S while their integrated forms are X®, X® and S respectively. denotes an n-particle 
tree-level colour ordered amplitude. 



with, 

/ d^/ifo = / ^LO + / ^NNLO (2-15) 

On the other hand, the double unresolved antenna functions X® in da^ b NLO and the 
colour-unconnected X®X® terms in da^^ LO can immediately be integrated over the phase 
space of both unresolved particles and appear directly in the m-particle final state, 

/ da NNLO = / d *NNLO + / d °NNLO- ( 2 ' 16 ) 
J2 J2 J2 

These integrated contributions will be discussed elsewhere. 

We now turn to a detailed discussion of each of the terms in Eq. (2.14). 

2.1.1 Cancellation of explicit infrared poles in da^j^ LO : Ji^nnlo 

In the antenna subtraction approach, the single unresolved configuration coming from 
the tree-level process with two additional particles, i.e. the double-real process involving 
(m + 4)-partons, is subtracted using a three-particle antenna function - two hard radiator 
partons emitting one unresolved parton. Once the unresolved phase space is integrated 
over, one recovers an (m + 3)-parton contribution that precisely cancels the explicit pole 
structure in the virtual one-loop (m + 3)-parton matrix element. 

The integrated subtraction term formally written as d&fj^ LO is split into three 
different contributions, depending on whether the hard radiators are in the initial or final 
state. 

When both hard radiators i and k are in the final state, then is a final-final 
(FF) antenna function that describes all singular configurations (for this colour-ordered 
amplitude) where parton j is unresolved. The subtraction term for this single unresolved 
configuration, summing over all possible positions of the unresolved parton, reads, 



d& NNLO = ^NNLO d $m+2(P3, • • • , Pm+A', Pi, P2) TT~ 

' Jm+2 

£ Xf jk \M m+ s(. . . , I, K, . . .)\ 2 Jt +1 \{pU +3 ). (2.17) 



perms 

X 

3 
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Besides the three parton antenna function X^- k which depends only on pi, pj and p^, the 
subtraction term involves an (m + 3)-parton amplitude depending on the redefined on- 
shell momenta pj and px, whose definition in terms of the original momenta are given in 
Appendix A.l. The (m+3)-parton amplitude also depends on the other final state momenta 
which, in the final-final map, are not redefined and on the two initial state momenta p\ 
and p 2 . This dependence is manifest as ellipsis in (2.17). The jet function is applied to the 
(m + l)-momenta that remain after the mapping, i.e. p%, . . . ,pi,px, ■ ■ ■ ,Pm+4- 

To perform the integration of the subtraction term in Eq. (2.17) and make its infrared 
poles explicit, we exploit the following factorisation of the phase space, 

d3> m+2 (f>3, • • • ,Pm+4;Pl,P2) 
= d$ m+ i(p 3 , . . . ,PI,PK, ■ ■ ■ ,Pm+4-,Pl,P2) ■ d$X ijk {Pi,Pj,Pk-,Pl + Pk) (2.18) 

dxi dx 2 

= d$ m+ i(p 3 , . . . ,p m+ 3] X 1 p 1 ,X 2 P2) 

Xl X2 

xS(l - xi)<5(l - x 2 )d$ Xijk (pi,Pj,Pk;pi +Pk), (2-19) 

where we have simply relabeled the final-state momenta in passing from (2.18) to (2.19). 
In (2.19) the antenna phase space d<3?x ijfe is proportional to the three-particle phase space 
relevant to a 1 — >■ 3 decay, and one can define the integrated final-final antenna by 

X? jk {s IK , Xl ,x 2 ) = ^ J d$ Xijk X% k 8{l - Xl )6(l - x 2 ) , (2.20) 

where C(e) is given in (2.10). 

The integrated single-unresolved contribution necessary to cancel the explicit infrared 
poles from virtual contributions in this final-final configuration then reads, 

f ,.S,a,(FF) K rR V ST [a* ( \ 1 ^Xidx 2 

/ da NNLO = N NNLO / d ®m+l{P3, ■ ■ ■ ,Pm+3! XlPl, X 2 p 2 )~ 

J I perms ^ ^ m+1 Xl %2 

x X U s ik^u^)\M m+3 (...,i,k,...)\ 2 Jt +1) ({pU +3 ) , (2.21) 

ik 

where the sum runs over all colour-connected pairs of final state momenta (pi,Pk) an d 
the final state momenta /, K have been relabelled as i, k. Expressions for the integrated 
final-final three-parton antennae are available in Ref. [21]. 

When only one of the hard radiator partons is in the initial state, Xij^ is a initial-final 
(IF) antenna function that describes all singular configurations (for this colour-ordered 
amplitude) where parton j is unresolved between the initial state parton denoted by i 
(where i = 1 or i = 2) and the final state parton k. The antenna only depends on 
these three parton momenta Pi,Pj and pk- The subtraction term for this single unresolved 
configuration, summing over all possible positions of the unresolved parton, reads, 

d& NNLO = ^NNLO d ^m+2(P3, ■ ■ ■ , Pm+4', Pi , P2) 

perms ° m+2 

x E Yx° jk \M m+3 (...J,K,...)\ 2 jt +1 \{ P } m +3). 

i=l,2 j 
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(2.22) 



As in the previous final-final case, the reduced (m + 3)-parton matrix element squared 
involves the mapped momenta I and K which are defined in Appendix A. 2. Likewise, the 
jet algorithm acts on the (m + l)-final state momenta that remain after the mapping has 
been applied. 

In this case, the phase space in (2.22) can be factorised into the convolution of an 
(m + l)-particle phase space, involving only the redefined momenta, with a 2-particle 
phase space [71]. For the special case i = 1, it reads, 

d$m+2(P3, • • • ,Pm+4-,Pl,P2) 

dxi Q 2 

= d$ m+ i(p 3 , . . . ,p K , . . . ,p m+ 4;x 1 p 1 ,p 2 ) S(xi - xi) -z-d$ 2 (Pj,Pk;Pi,q) (2-23) 

= d$ m+ i(p 3 , . . . ,p m+ 3;x 1 p 1 ,x 2 p2) — — S(xi - xi)S(l - x 2 )Q-d$2(pj,Pk;pi,q) , 

X\ X 2 27T 

(2.24) 

with Q 2 = —q 2 and q = pj + p^ — p\ and where we have also relabelled the final-state 
momenta in passing from (2.23) to (2.24). The quantity x\ is defined in Eq. (A. 7). 

Using this factorisation property one can carry out the integration over the unresolved 
phase space of the antenna function in (2.22) analytically. We define the integrated initial- 
final antenna function by, 

X? dk {s lK ,x 1 ,x 2 ) = -^J d^S(x 1 -x 1 )S(l-x 2 )X° lijk , (2.25) 

where C(e) is given in (2.10). Similar expressions are obtained when i = 2 via exchange of 
x\ and X2, 

X2jk(*2K,*i,X2) = ^ J d^S(x 2 -x 2 )S(l- Xl )X^ jk . (2.26) 

The explicit poles present in da§^ LO (defined in (2.13)) associated with the colour- 
connected initial-final radiators i and k can therefore be removed with the following form, 



f ,-S,a,(IF) » rRV f A ih I \ 1 dx\dx 2 

da NNLO =Nnnlo >^ d$ m+l (p 3 ,. . . ,p m+3 ;xipi,x 2 p2)- d 

h p j^ s J S m+1 Xi X 2 

x E E ^ fc (^^i^2)|^ m+ 3(...,/,fc,...)| 2 Ji m+1) (M m+3 ). (2.27) 

i=l,2 k 

In this expression, only the redefined momentum K has been relabelled (K — > k). The 
rescaled initial state radiator /(= i) is not relabelled and appears in the functional depen- 
dence of the integrated antenna and in the matrix-element squared. Explicit forms for the 
integrated initial-final three-parton antennae are available in Ref. [71]. 

If we consider the case where the two hard radiator partons i and k are in the initial 
state finally, then X^j is a initial-initial (II) antenna function that describes all singular 
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configurations (for this colour-ordered amplitude) where parton j is unresolved. The sub- 
traction term for this single unresolved configuration, summing over all possible positions 
of the unresolved parton, reads, 

d °NNLO = NnNLO E d$ m+2 (p3, . . . ,p m+ 4-,Pl,P2) 

Orn-i-2 

perms i>»-r^ 

x E T, X kj\ M ^(...J,K,---)\ 2 Jt +1) (P3,...,P m+ 4). (2.28) 

i,fc=l,2 j 

where as usual we denote momenta in the initial state with a hat. The radiators i and k 
are replaced by new rescaled initial state partons I and K and all other spectator momenta 
are Lorentz boosted to preserve momentum conservation as described in Appendix A. 3. 

For the initial-initial configuration the phase space in (2.28) factorises into the con- 
volution of an (m + l)-particle phase space, involving only redefined momenta, with the 
phase space of parton j [71] so that when i = 1 and k = 2, 

d$ m+2 (f>3, • • • ,Pm+4;Pl,P2) 

dx\ dx2 



d$ m+ i(p 3 , . . . ,p m +4;xipi,x 2 p2) x X\X 2 5{x\ - X\) 5(x 2 - x 2 ) [dpj] 

X\ X2 



(2.29) 



where the single particle phase space measure is [dpj] = d d ~ 1 pj /2Ej /(27r) d ^ 1 and xi is 
defined in Eq. (A. 9). 

The only dependence on the original momenta lies in the antenna function j and 
the antenna phase space. One can therefore carry out the integration over the unresolved 
phase space analytically, to find the integrated antenna function, 

x i2,j( s l2> x u x 2) = -^rr J[dp j ]x 1 x 2 d(x 1 -x 1 )S(x2-X2)X^ j , (2.30) 

where C(e) is given in Eq. (2.10). Explicit forms for the integrated initial-initial three- 
parton antennae are available in Ref. [71]. 

We can therefore remove the explicit poles present in da^ LO (defined in (2.13)) 
associated with the colour-connected initial-state pair i and k with the subtraction term, 



I d °NNLO = Nnnlo yZ I d® m+l (p 3 , . . . ,p m+3 ;x 1 pi,x 2 p2)-^- 



dxi dx2 



|>rni]s . Wl x l x 2 

X E </^,^,^)l^ m +3(...,/;^,...)| 2 4 m+1) (Mrn +3 ). 
i,k=l,2 

(2-31) 

As in the initial-final case, the redefined initial state momenta I = i and K = k are not 
relabelled, neither in the functional dependence of the integrated antenna Xf k Xj, Xk) 
nor for the reduced matrix element squared |A4 m +3(. . . I, K, . . .)| 2 . 

Each of the integrated antenna have an explicit dependence on the variables x\ and 
x 2 as stated in eqs. (2.20), (2.25) and (2.30). 
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Summing over the different final-final, initial-final and initial-initial configurations, we 
find that the explicit poles in da§^ LO are removed by j l da S ^ a NLO to yield an integrand free 
from explicit infrared poles over the whole region of integration. This is of course merely a 
consequence of the cancellation of infrared poles between a virtual contribution and inte- 
grated subtraction terms from real emission. However, as discussed later in Section 2.2.2, 
d&N^NLO an d Ji d<7jV7VLO develop further infrared singularities in singly unresolved regions 
of the phase-space which do not coincide. Therefore, we have to introduce a further sub- 
traction term, da^ s NLO , to ensure an integrand that has no explicit global e-poles, and that 
does not have implicit singularities in single unresolved regions. 

2.1.2 Integration of iterated antennae: f t da^Nix) 

As discussed earlier in Sect. 2.1, contributions to the double real subtraction term due to 
colour-connected or almost colour-unconnected hard radiators that have the generic form 
X 3 x X 3 must also be added back integrated over the phase space of the first (outer) antenna. 
The structure of these terms is very similar to the single unresolved contributions of the 
previous sub-section; each term with an "outer" final-final antenna present in da^^o 
produces a contribution given by, 

1 dxi dx 2 



J i da NNLO FF) = NnNLO j d<f> m +l(P3, ■ ■ ■ ,Pm+3;X 1 p 1 ,X 2 P2) 



Sm+l Xi X 2 



x X? jk (s tk ,x 1 ,x 2 )Xl({p} m+3 )\M m+2 ({p} m+ 2)\ 2 4 m) (M m+2 ) • 

(2.32) 

In this case, i and k represent the momenta in the set {p} m +3- The set of momenta denoted 
by {p}m+2 is obtained from {p} m+3 set through a phase space mapping that is determined 
by the type of the unintegrated "inner" antenna X® which is fixed by the corresponding 
term in da% NLQ . 

When the outer antenna is of initial-final or initial-initial type, we find the integrated 
forms, 



J da NNLO IF) = ^NNLO j d$ m+1 (p 3 , . . . ,p m+3 ; X X p U X 2 p 2 ) 



1 dx\ dx 2 



S m +i x\ x 2 

X ^ifc(Sifc^l^2)^({p} m+3 )|A4 m+2 (M m +2)| 2 «/i m) (Mm+2) , 

(2.33) 

[ ,~S,(b,c),(ii) Kr Rv f ,^ 1 \ f dxi dx 2 

da NNLO = N NNLO / d$ m+ i(p 3 , • • • ,Pm+3; X\p\, X 2 p 2 )- 

Jl J dm+1 X\ X 2 

x ^ J (^,x 1 ,x 2 )X 3 (M m+3 )|A4 m+2 (M m+2 )| 2 J^{{p} m+2 ) • 

(2.34) 

S h 

Note that the (b, c)-type terms reflect different physical origins; da^ NLO is designed 
to account for the unresolved contributions from two colour-connected unresolved par- 
tons, while do S pf NLO treats the singularities from two almost colour-connected unresolved 
partons. 
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S b 

The general structure of da^ NLO consists of groups of terms involving the difference be- 
tween a four-parton antenna and products of three-parton antennae. The latter removes the 
single unresolved limits present in the former. There are two distinct types of four-parton 
antenna which, for sake of argument, we label X®(hi, si, 82, ^2) and X®(hi, s±, S2, ^2)) where 
hi, ti2 represent the hard radiators and si, S2 the unresolved particles. By construction, 
there are no unresolved limits when h\ or h<i become unresolved. In both, the four partons 
are colour-connected. The key difference between X\ and X\ is that in X®, si (52) is 
colour-connected to h\ (J12) and S2 (si), while in X®, si (^2) is colour-connected to hi and 
h 2 . 

In the case of X®, da^J b NLO has the form, 
X®(hi, si, s 2 , h 2 )\M n +2(- ■ ■ , (hisis 2 ), (siS2h 2 ), ■ ■ -)\ 2 

-Xl(h U S U S 2 )X°({h^), {^2)M)\M n (. . . , ((/il^XTlSa)), ((sTSM, • • -)| 2 

-Xl( Sl , s 2 , h 2 )Xl(h 1 , (^S), (S^))\M n (. . . , (M^7 2 )), • • -)| 2 (2-35) 

where the two X® x A3 terms subtract the single unresolved limits present in X®. Most 
antennae fit into this class including, for example, the quark-antiquark antenna A®, the 
quark-gluon subantenna Z?4 a [58] and the gluon-gluon subantenna F® a [75]. Upon integra- 
tion, the A3 x A3 terms produce integrated antenna that cancel against the explicit poles 
present in A3 that will be discussed in Section 2.2.2. 

Similarly, one must also subtract the limits where S\ and S2 are unresolved from X®, 

X2(hi,si, s 2 , h 2 )\M n+ 2{- ■ • , hisis 2 , s 1 s 2 h 2 , • • -)| 2 
-Xl(h u s 1 ,h 2 )X$((h^i), 8 2 , (sM)\M n {. ■ ■ , (h^fs 2 ), (s^M), . . .)| 2 
-Xi(h 1 ,s 2 , h 2 )X$((h^), si, (SS))|M n (. . • , (h^fs!), (s^2~h 2 )), ■ ■ -)| 2 - (2-36) 

Antennae of this type include the quark-antiquark antenna A®, the quark-gluon subantenna 
D\ b [58] and the gluon-gluon subantenna F® b [75]. As before, after integration the A3 x A3 
terms produce integrated antenna that cancel against the explicit poles present in X\. 
However, unless the single unresolved contribution d<r^ LO has the form, 4 

X%(h u si, h 2 )\M n+ i(. . . , (mT), 8 2 , SS), . . .)| 2 , (2.37) 

it is clear that Eq. (2.36) oversubtracts the double unresolved limits. In this 
correction term of the form 

+^(/ ll)S i,/ l2 )i 3 o (M ) 82, Smim„(. . . , ((mI> 2 ), {S2IM2)), ■ ■ 01 2 

+±Xl(h 1 ,S2, h 2 )Xl{{h^2), Sl , . . , ((Qsi), (s^h 2 )), ■ ■ -)| 2 , (2-38) 

must be added to (2.36) to compensate. In the case of A4, da^ b NLO is obtained by summing 
Eqs. (2.36) and (2.38). The terms in (2.38) are identified as having repeated hard radiators 



4 This typically happens when both si and S2 are photon- like, see for example Eq. (8.2) of Ref. [58]. 
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i.e. the hard radiators of the outer antenna (hi and h 2 ) produce mapped momenta (his) 
and (sh 2 ) that become the hard radiators of the inner antenna. They can be integrated 
according to Eqs. (2.32), (2.33) and (2.33), to give contributions of the form, 

+^(s hlh2 ,x 1 ,x 2 )Xi(h 1 ,s 2 , h 2 )\M n {. . . , (SS), (ST), . . .)| 2 . (2.39) 

These terms produce poles that cancel against the integrated wide angle soft terms dis- 
cussed in Section 2.1.3. and the real-virtual subtraction terms discussed in Section 2.2.3. 
The presence of D® b and F® b antennae in dafj NLO is an indicator that there are contribu- 
tions from wide angle soft radiation, see for example Refs. [58,75]. 

Similarly, the almost colour-connected contributions in da^^LO are produced by ma- 
trix elements of the type \M n + 2 (. ■ ■ , a, si, hi, s 2 , h 2 , . . .)| 2 and have the form [21,75], 

--Xl(a, si,hi)Xl((hiSi), s 2 , h 2 )\M n {. . . , (asi), {(h lSl )s 2 ), (s 2 h 2 ), . . .)| 2 

~Xl{hi,s 2 , h 2 )Xl{a, si, (Ms7))\M n (. ...(osT), ((£^1), (ST), . . - )| 2 - (2-40) 

Note that unlike the contributions described in Eq. (2.38), the radiators for the inner 
antenna are not both constructed from mapped momenta. Upon integration, these terms 
also produce poles that cancel against wide angle soft terms and will be further discussed 
in Section 2.2.3. 

2.1.3 Integration of large angle soft terms: J x da^^LO 

For processes involving soft gluons the real-real channel has an additional subtraction 
contribution denoted by daff e NLO due to large angle soft radiation [58,59,62,75]. This 
term removes the remnant soft gluon behaviour associated with the phase space mappings 
of the iterated structures of the double-real subtraction contributions da^^'u)- ^ n these 
contributions two successive mappings have been applied and the gluon emission can occur 
before or after the first mapping. The large angle soft subtraction terms is constructed 
to account for that and the soft eikonal factors can involve either unmapped or mapped 
momenta. 

The wide angle soft terms are however not uniquely defined, only their soft behaviour 
is determined. They can be obtained after two successive mappings of the same kind have 
been applied. In other words, either two final-final, two initial-final or two initial-initial 
mappings are used. Those are presented in [75]. 

Alternatively, they can be constructed after the application of two successive mappings 
where the first of those is always a final-final mapping, (i,j,k) — > (I,K). To achieve this 
we modified the large angle soft contribution of the double-real piece, given by Eqs. (5.9), 
(5.13) and (5.17) of Ref. [75], to the equally valid forms given in (C.l), (C.2) and (C.3) 
respectively. Following this latter choice, the integrated large angle soft subtraction terms 
can simply be obtained through integration over the final-final antenna phase space d&x ijk 
given in eq. (2.19). 

The unintegrated soft factor terms are eikonal factors of the form 

q 2s ac 
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where a, c denotes initial or final state external partons. The hard radiators, a and c, may 
be momenta present after the first mapping, or momenta produced by the second mapping, 
i.e. 

Pa, Pc = Pa, Pc e {p} m+ 3 or p a , p c = p a < , p c > G {p} m +2- (2.41) 

Because j is eliminated in the first mapping, it is not a member of {p} m + 3 or {p} m +2- 
Nevertheless, in the latter case, the second mapping should also be applied to pj to produce 
a new soft momentum py . When this is a final-final or initial-final mapping, it is a trivial 
mapping pj — > p,y = pj. However, when the second mapping is of the initial- initial form, 
then pj —7- pji = pj as given by Eq. (A. 8). 

In general, and after factorisation of the phase space according to Eq. (2.18), one 
always finds differences of terms of the form, 

rRR *~ • 



N NNLO d®x ijk d$ m+ i(p 3 , . . . ,p m+3 ;p 1 ,p 2 ) 



>Jm+l 

X \ [Sajc ~ S alj , d ] X$({ P } m+3 ) \M m+ 2({p} m+ 2)\ 2 jt ] ({p}m + 2) , (2.42) 

where X® is the antenna corresponding to the second mapping, i.e. final-final, initial-final 
or initial-initial. 

For the final-final kinematical configuration, where the second mapping transforms the 
momenta (I, I, K) — > (/' ,K'), we find a contribution of the form, 

/ &° S NNLO ] '=^NNLO [ ^X ijk d$ m+1 (p 3 , . . . , Pm+3 ; Pl , p 2 ) — 



X - S l'iK' ~ S IjK ~ S ajl' + S ajl ~ S K 'jb + S Kjb 



^IlK 



1 

2 Y^I'jK' ~ ^ l j J^. - >~> a jf t- ^ajl - ^K'jb 

,Pji ,p K i , . . . ,p m+3 ) . (2.43) 



We see that there are three pairs of soft antennae which correspond to inserting the soft 
gluon around the radiator gluons, in this case, between a and I', between I' and K' and 
between K' and b. 

Similarly, when the second mapping is of the initial-final type corresponding to the 
mapping (n,l,K) — > (N,K'), we have, 



/ d °NNLO =^NNLO [ d ®X ljk V d$ m+i (p 3 , . . . , p m+3 ; Pl , p 2 ) 
J J m+l bm + 1 

X ^ o [^JVj-K"' ~~ ^ K ~ ^K'jb + SKjb - S a + S a jh X® lK 
n=l,2 

x \M m+2 (. ..,a,N,K',b,...)\ 2 (p 3 , ■ ■ ■ ,Pk>,- ■ ■ ,P m+ z) • (2.44) 



Once again, there are three pairs of soft terms, corresponding to inserting the soft gluon 

in the three positions into the colour-connected gluons, . . . , a, N, K' , b, In this case, 

the last two terms simply cancel. This is because the radiator momenta appears in both 
the numerator and denominator of the eikonal factor, coupled with the fact that the initial 
momentum is simply scaled in the initial-final mapping, S a jn = S a j^. 
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Finally, when the second mapping is of the initial-initial type: (n,p, /, . . . , a, b, . . .) 
(N, P, . . . , a, b, . . .) we have, 

/ d& NNLO = Nnnlo \ d$x yfc V d$ m+ i(p3,...,p m+3 ;pi,p2)-^— 
7 7 m+l 



X 



NjP 



Pjb 



X 



np,i 



n,p=l,2 

x |A^ m+2 (. . . , 5, JV.P, 6, . . .)| 2 4 m) fe, • • • ,Pm +3 ) 



(2.45) 



As anticipated, when the soft factor involves hard radiators produced by the initial-initial 
mapping, the soft momentum pj is replaced by the momentum p~- obtained by applying 
the boost for the initial-initial mapping directly to momentum pj. 

When the unresolved momentum is unboosted and denoted by pj , the integrated form 
of the soft factor is given by, 



S(Sac,SlK,Xac,IK) = J d®X ijk S a jc, 



(2.46) 



where a and c are arbitrary hard radiator partons. After relabelling the momenta such 
that I,K ->■ i, k, then [58], 



where we have defined 



1 1 



ln(xac,ik) - Li 2 - 



1 — X 



ac,ik 



%ac,ik 



111 

12 



+ 0(e) 



Sac^ik 



(2.47) 



(2.48) 



Since the soft eikonal factor is invariant under crossing of one or two partons from the final 
to the initial state the integrated soft factor is independent of where the external partons 
are situated. 

The integrated final-final soft factor involving a boosted unresolved momentum (de- 
noted by p~- in the soft eikonal factor is defined by, 

1 



>5(s ac , Sik, x aCi jfe) — 



C{e) 



d ®x lJk S a]c . 



(2.49) 



Since S a bc is composed of Lorentz invariants, we can simply invert the Lorentz boost that 
mapped j — > j (A. 8) such that 

S a]c = Sajc (2.50) 
where a, c are also obtained by inverting the same Lorentz boost so that 

1 



C(e) 



d$X ijk Sajc- 



(2.51) 



The RHS has precisely the same form as Eq. (2.46) with p a — > p^ and p c — > pc, so that 

•S{s ac , Sik, x ac ^k) = S{s ac , Sjk, x ac jk). (2.52) 
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Furthermore, we can exploit the fact that S only depends on invariants and can apply the 
same Lorentz boost such that 



so that, 



{Pa,Pc,Pi,Pk} -> {Pa,Pc,Pi,Pk} 



'Si.Saa $iki X a c,ik) — S(^S ac , S^, X QC 



(2.53) 



The integrated soft factors only contribute to the soft region (xi = x 2 = 1)- Using 
the (to + 1) phase space factorisation as given in Eq. (2.18) and inserting the integrated 
large angle subtraction factors into Eqs. (2.43), (2.44) and (2.45) while relabelling I' — > I, 
I -> i, K' -> K and K -)• k we find, 



/ ^NN^P = ^NNLO E / d$m+l(P3, • • • ,Pm+^X 1 p 1 ,X 2 P2, 

m+1 



1 dxi dx2 



Sm+l %1 X 2 



x-Xf lk 5(l- Xl )S(l-x 2 ) 



S(sjk, Sik, XlK,ik) — 5(Sjfc, Sit, Xik,ik) ~ 5(s a /, Si k ,X a l 



ik) 



+S{s a i, Sik,X ai)i k) — S(sxb, Sik,XKb,ik) + 5(Sfcfe, S 

ik} Xkb^ikj 

x\M m+2 {. ..,a,I,K,b,...)\ 2 4r } (p 3 , ...,PI,PK,-.. ,Pm+3) , 

/..Se(iJf) . .pi/ v-^ f / \ 1 dxi dX2 
da NNLO = N NNLO / d$ m+ i(p 3 , . . . ,p m+3 ; XiPi, X 2 p 2 ) ~ 

m+1 J m+1 Xl X2 

X E " £!H(1 - X 2 ) 



(2.54) 



n=l,2 



^^■ S NK'> ^iki ^NK,ik) ^i.^nki Siki Xfik,ik) 



-S(s Kb , s ik ,x K b,ik) + 5(s fcb , s ik ,x kbyik ) 
x |-M m+2 (- ■ ■ ,a,N,K,b, . . .)| 2 Ji m) (Ps, ...,Pk,--- ,Pm+ 3 ) , 



(2.55) 



, ~S,e,(//) _ /w-W 
NNLO ~ yV NNLO 



E / d$ m+ i(p 3 , . . . ,p m +3;XlPl,X 2 P2) 
m+1 ^ 



1 dxi dx 2 



E ^,*(l-xi)«(l-x 2 ) 



n,p=l,2 



5 m+ i Xl x 2 
_ ^(^JV\P' ~^NP~ik) + ^(^npi ^ik-i Xfip^ik) 



+'"' ( S aA r ' ^aW" ik) ^{ s ahi s iki X a h,ik) 4" ^{ s pbi s i k i x pb ifc) S{ s pbt s iki Xpb,ik) 

x\M m+2 {. ..,a,N,P,b,.. .)| 2 jt ] (P3, ■ ■ ■ ,j5m+ 3 ) • (2.56) 

Note that the double e-poles cancel within the square brackets, so that, for example, the 
leading pole for the square bracket in Eq. (2.54) is 

SaiSlKSkb\ 



1 



log 



+ 0(1) 



2e \SalSi k SKbJ 

and similarly for Eqs. (2.55) and (2.56). 

Moreover, this combination of integrated soft factors does not have any single unre- 
solved limit and therefore does not require further subtraction. 
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a a b,c 


drr VS 
NNLO 

h NNLO 


n\M° m+2 \ 2 Xl\Ml +2 \ 2 X°X$\M° m+2 \ 2 
^\Ml +2 ? Xl\M] n+2 ? X§X§\Ml+ 2 \* 



Table 2: Type of contribution to the real- virtual subtraction term d<7^f rLO , together with the 
integrated form of each term. The unintcgratcd antenna functions are denoted as X® and X% while 
their integrated forms are X§ and X\ respectively. I-M* | 2 denotes the interference of the tree-level 
and one-loop n-particle colour ordered amplitude while |.M°| 2 denotes the square of an n-particle 
tree-level colour ordered amplitude. 



2.2 Subtraction terms for one-loop single-unresolved contributions: do^ s NLO 

As we have seen before, the contribution obtained by integrating the single unresolved 
subtraction terms present in d<7^£ precisely cancels the pole structure of the one- 
loop (m + l)-parton contribution daff^ L0 . The real- virtual subtraction term da^ s NL Q 
of Eq. (2.12) therefore has to subtract three types of contributions: 

(a) Single unresolved limits of the virtual one-loop (m + 3)-parton matrix element. 
Each subtraction term takes the form of a tree- or one-loop-antenna function multi- 
plied by the one-loop or tree- colour ordered matrix elements respectively which we 
denote by dd^^ a L0 . 

(b) Terms of the type X® X® that cancel the explicit poles introduced by one-loop matrix 

elements and one-loop antenna functions present in da^f^ a L0 . This term is named 
, - VS,b 
® NNLO' 

(c) Terms of the type X® X® that compensate for any remaining poles, specifically those 
produced by J da^^'ib anc ^ / ^nnlo-> an< ^ wmcri are denoted by da^f^Q, 



so that 



A-VS _ A -VS,a . ,~VS,b , A *VS,c / K7 \ 

aa NNLO — aa NNLO + aa NNLO + aa NNLO' l z - D 'J 



The types of contributions present in each of these terms are summarised in Table 2. 

Note that da\j S NLO is a subtraction term, and it must therefore be added back in 
integrated form to the m-parton final state. There it combines with the twice integrated 
double real subtraction term J 2 da^ NLO and the mass-factorisation counterterm denoted 
by da^KLo according to Eq. (1.5) to provide the singularity structure necessary to cancel 
the pole structure of the double virtual contribution da)^ L0 . 

2.2.1 One-loop single-unresolved contributions: da^^LO 

In single unresolved limits, the behaviour of (m+3)-parton one-loop amplitudes is described 
by the sum of two different contributions [78-82]: a single unresolved tree- level factor 
times a (m + 2)-parton one-loop amplitude and a single unresolved one-loop factor times a 
(m+2)-parton tree- level amplitude, as illustrated in Figure 1. Accordingly, we construct the 
one-loop single unresolved subtraction term from products of tree- and one-loop antenna 
functions with one-loop and tree-amplitudes respectively. 
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Figure 1: Illustration of NNLO antenna factorisation representing the factorisation of both the 
one- loop "squared" matrix elements (represented by the white blob) and the (m+ l)-particle phase 
space when the unresolved particles are colour-connected. The terms in square brackets represent 
both the three-particle tree-level antenna function X^- k and the three-particle one-loop antenna 
function X\- k and the antenna phase space. 
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j 



Analogously to the case of handling explicit infrared poles present in J x da^-^LO or 
to deal with single unresolved limits of tree level amplitudes, we need to decompose this 
contribution into three parts according to where the two hard radiators are situated. We 
note that the momentum mappings which implement momentum conservation away from 
the single unresolved limits of each configuration are the same as in the previous section, 
namely equations (A. 4), (A. 6) and (A. 8). 

In the final-final configuration, the subtraction term is given by, 

,,VS,a,(FF) KrRV ST^ [ax, ( \ 1 dx\dx 2 

X? jk 8(1 - Xl )S(l - x 2 ) \Ml n+2 (. ..,I,K,...)\ 2 J^({p} m+2 ) 

+ X} jk 8(1 - Xl )8(l - x 2 ) \M m+2 (. ..,I,K,...)\ 2 jH({p} ro+2 ) 

(2.58) 

where the one-loop three-parton antenna function Xj- k depends only on the antenna mo- 
menta Pi,Pj,Pk- X-ijk correc tly describes all simple unresolved limits of the difference 
between an (m + 2)-parton one-loop corrected squared matrix element and the product 
of a tree-level antenna function with the m-parton one-loop corrected squared matrix el- 
ement. It can therefore be constructed out of colour ordered and renormalised one- loop 
three-parton and two-parton matrix elements as 

x m ~ s w ^rj 2 ~ x ijk , (2.59) 

where SijkjK denotes the symmetry factor associated with the antenna, which accounts 
both for potential identical particle symmetries and for the presence of more than one 
antenna in the basic two-parton process. 

It should be noted that X}-, is renormalised at a scale corresponding to the invariant 
mass of the antenna partons, Sij k , while the one-loop (m + 2)-parton matrix element is 
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renormalised at a scale /i 2 . To ensure correct subtraction of terms arising from renormali- 
sation, we have to substitute 



in (2.58). The terms arising from this substitution will in general be kept apart in the 
construction of the colour ordered subtraction terms, since they all share a common colour 
structure /3q. 

Similar subtraction terms are appropriate in the initial-final and initial-initial config- 
urations. In the initial-final case, we have, 



d ^NNLO F) = J^NNLO E / d< ^m+l(P3, ■ ■ ■ ,Pm+3; X 1 p 1 ,X 2 p 2 ) 



1 dxi dx 2 



S m+ i xi x 2 



* EE 

i=l,2 j 



perms 

X° jk 5(1 - Xl )S(l - x 2 ) \Ml +2 {. . . , I, K, . . .)\ 2 JH ( {p}m+2 ) 

+ Xl Jk 5(1 - Xl )5(l - x 2 ) \M m+2 (. ..J,K,...)\ 2 Jt ] (Mm+2) 

(2.61) 

where the one-loop antenna function X\- k is obtained from Xj- k by crossing parton i from 
the final state to the initial state. 

The initial-initial subtraction term is 

d(J NNLO = N NNLO I. I d$ m+l(P3, • • • ,Pm+^ XlPl, X 2 p 2 )- 

„t? ns J b rn+l Xi X 2 

X&j 5(1 - Xl )5(l - x 2 ) \Ml +2 (. . . , /, K, . . .)| 2 Jt\{p} m+ 2) 



x E E 

i,k=l,2 j 



+ X} Kj 5(1 - Xl )5(l - x 2 ) \M m+2 (. ..J,K,...)\ 2 jH({p} m+2 ) 

(2.62) 

with Xf k j obtained by crossing partons i and k from X\- k . We note that in order to 
fulfil overall momentum conservation, the initial-initial momentum mapping requires all 
the momenta in the arguments of the reduced matrix elements and the jet functions to be 
redefined. 

The one- loop single unresolved subtraction terms, X\- k , X\ - k and X\ k ■ can never be 
related to integrals of tree-level subtraction terms. Therefore, after integration over the 
three-parton antenna phase space, this component of the subtraction term must be added 
back in integrated form to the terms yielding the m-parton final state contribution. This 
can be accomplished using the techniques described in [91] and analytic expressions for all 
of the integrated three-parton one- loop antennae are available in Refs. [21,72,74]. 
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2.2.2 Cancellation of explicit infrared divergences in da^f^ a LO : d&^f LO 

By construction daj^^Lo correctly approximates the one-loop (m + 3)-parton contribution 
to (m + l)-jet final states in all single unresolved limits. However, the one-loop antenna 
functions and one-loop reduced matrix elements present in Eqs. (2.58), (2.61) and (2.62) in- 
troduce explicit e-poles which do not correspond to poles of the real-virtual matrix element 
(2.13) and therefore have to be removed. 

It is the purpose of this subsection to present a new term denoted by da^ff LO which 
removes the explicit e-poles present in dd^ S ^ a L0 . To achieve this, we introduce further 
partially integrated subtraction terms built with products of an unintegrated and an in- 
tegrated three parton antenna function times a reduced matrix element squared. For a 
given configuration, the integrated antenna will be of a given type (final-final, initial-final 
or initial-initial) multiplying a sum of three unintegrated antennae of each configuration 
type. 

Let us consider first the explicit poles in the one-loop reduced matrix elements present 
in Eqs. (2.58), (2.61) and (2.62). The pole structure of \M.] n+2 {{p} m +2)\ 2 is well understood 
as a sum of I*- 1 -* operators [76]. These e-poles can be simply subtracted using integrated X® 
antennae. To derive the relevant term, we simply add a contribution constructed from each 
1(1) 

operator present in the pole structure of \-Mj n+2 ( K {p} m +2)\ 2 according to the relation 
shown in Table 3, 

2I^(s ab )^ 1 Lxl (2.63) 

^ab 

The particle type (x, y) fixes the flavour of antenna, while the momenta (a, b) determines 
whether the integrated antenna is final-final, initial-final or initial-initial. This is sum- 
marised by adding the contribution obtained by the replacement 

I^+2({PW2)I^(1 " Zl)<5(l - X 2 ) -> (^x^(s ab , Xl ,X2)) \M° m+2 ({p} m+2 )\ 2 , 

ab ^ab / 

(2.64) 

in Eqs. (2.58), (2.61) and (2.62). Here ab runs over the pairs of colour-connected particles, 
and X® (s a b) is an integrated antenna that depends on the particle type and on whether 
a and b are in the initial or final states. This contribution guarantees that the e-poles 
associated with soft and final state collinear singularities introduced through the one-loop 
matrix elements in Eqs. (2.58), (2.61) and (2.62) are cancelled. 5 Note that the momenta 
a and b lie in the set {p} m +2 and that the invariant mass associated with each antenna is 
constructed from mapped momenta and cannot be related to terms coming from J l da 5 ' 1 . 

There are also explicit poles present in the one-loop antenna, Xj- k which are produced 
by the physical matrix elements making up the antenna (2.59). These can also be described 
by integrated antenna, however now the relevant mass scale is constructed from momenta 
before the mapping i.e. the momenta lying in the set {p} m +3- 

The first type of poles are associated with l-M^J 2 in Eq. (2.59). The second contribu- 
tion comes from the term proportional to l-Mj^l 2 . For example, the explicit poles in Xj- k 



3 Any remaining initial state collinear singularities are removed by the mass factorisation counterterm. 
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Table 3: The relation between iW operators and integrated antenna showing the dependence on 
whether the pair is in the final-final, initial-final or initial-initial state. 
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Table 4: Number of colour-connected pairs Nx in the one- loop antenna X\, and the coefficient 
M x . 



are removed by adding the contribution obtained by the replacement 
X} jk 5{\-x 1 )5{\-x 2 )^ 

^2 -^x x i( s cd,xi,x 2 ) J - ^-X$(s IK ,xi,x 2 



cd cd 



ik 



ijki 



(2.65) 



where cd runs over the Nx pairs of colour-connected particles within the antenna. Nx 
and Mx depend on the type of antenna, and their values are given in Table 4. Similar 
expressions are valid for initial-final and initial-initial antenna. The terms in the sum in 
(2.65) are systematically provided by daffff^Q, *' e " ^ ne -^3 x -^3 type terms coming 
from Eqs. (2.35) and (2.36). This is another example of the cancellation of infrared poles 
between a virtual contribution and integrated subtraction terms from real emission. On 
the other hand, the last term must be introduced as a new subtraction term in da^^f L0 . 



V S b 

To be more explicit, in the final-final configuration, da NN . LO is given by, 



da 



VS,b,(FF) 
NNLO 



K rRV 
JV NNLO 



E 

perms 



d$ 



m+l(P3 



,p m+ 3;x 1 p 1 ,x 2 p2) 



dx± dx 2 



Wl x l x 2 



E 



X] 4x X 3( S ab,Xl,X 2 ) j - ^^Xi(s IK ,Xi,X 2 ) 
V ab b ab J b ik 

x\M° m+2 ({p} m+2 )\ 2 JtHMm+2). 



X 



ijk 



(2.66) 



Here pi and px are the momenta produced by the mapping for the ijk antenna, while ab are 
the pairs of colour-connected particles appearing in the matrix element A^^ +2 ({p} m +2)- 
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Similar subtraction terms are appropriate in the initial-final and initial-initial configura- 
tions, 



da 



VS,b,{IF) 
NNLO 



NnNLO ^2 J d$ m+ i(p 3 , • • • ,Pm+3-,X 1 p 1 ,X 2 P2) 



1 dxi dx 2 



perms 



S m+ i Xi X 2 



i=l,2 j 



J2^^(s ab ,x 1 ,x 2 )) - —^-X$(s- iK ,x 1 ,x 2 ) 



ab ab / ik 

x|^ m+2 (M m+2 )| 2 Ji m) ({p} m+2 ) , (2.67) 

1 dxi dx 2 



^NNLO^ =NnNLO ^2 J d$ m+ i(p 3 , • • • ,Pm+Z-,XiPi,X 2 p 2 ) 



perms 



S m+ i xi x 2 



i,fc=l,2 j 



ik,j 



x\M m+2 ({p} m+2 )\ 2 Jt\{p} m+ 2) ■ 
2.2.3 Compensation terms for remaining poles: da^f^ L Q 



(2.68) 



In Sect. 2.1.2, we identified the integrated form of the double real radiation subtraction 
term that is due to almost colour-connected hard radiators. In the final-final case, this has 
the form, 

/ da NNLO = N NNLO / ^ m+1 {p 3 , . . . ,p m+3 ; ZlPl, X 2 p 2 )- 

Jl J &m+l Xi X 2 

x ^4 fc (^,x 1 ,x 2 )x 3 (M m+ 3)|>( m+2 (M m+2 )| 2 Ji m )(M m+2 ) 

(2.69) 

where the invariant mass is constructed from momenta in the set {p} m+3 . For each term 
like this, we introduce a subtraction term (with the opposite sign) where the appropriate 
invariant mass is constructed out of momenta in the set {p} m+2 i.e. the momenta produced 
by the mapping appropriate to the type of antenna X®({p} m+ 3). That is, 

,~VS,c,(FF) . rRV , 1 dx 1 dx 2 

da NNLO = -NnNLO / d$ m+ i(p 3 , . . . ,p m +r,XlPl,X 2 p 2 )- 

J im+1 X\ X 2 

x ^4 fc (s/x,x 1 ,x 2 )x 3 (M m+ 3)|^ m+2 (M m+2 )| 2 4 m )(M m+2 ) . 



Taken together, we find that 



f ,*S,c,(FF) . ,.VS,c,(FF) . r RV f ,^ / > 
J da NNLO + da NNLO = ^NNLO J d$ m+ l {P3, ■ ■ ■ , Pm+3] XlPl , X 2 p 2 , 



(2.70) 



1 dxi dx 2 



5 m+ i xi x 2 



Xijk( s ik,X!,X 2 ) - Xf jk (s IK ,xi,x 2 ) 



Xi({p}m +3 ) 



x \M m+2 ({ P } m+2 )\ 2 JtH{p} 



m+2) 



(2.71) 
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where the term in square brackets does not have a double pole in e, but has a leading 
singularity of the form 

Combinations of terms like this, plus similar contributions from the initial- final and initial- 
initial antenna, 

S,cXIF) ^ A ,VSMIF) _ K rRV__ [ Ath , „ ._ _J_ ^2 



[ ,~S,C,(IF) . ,.VS,C,{IF) .rRV ( Aih I ^ 1 

da NNLO + da NNLO = N NNLO d <$> m+ i{p 3 , . . . ,p rn+3 ;Xip 1 ,X 2 P2)-£- 
Jl J <->m- 



m+l X\ X2 



X 2 



^i,jk( s ik,xi,x 2 ) - X° jk (s- iK ,X!,X2) 



X° 3 ({p} m+3 ) 



x |X m+2 (M m+2 )| 2 4 m) (W m+2 ), (2.72) 

dxi dx2 



I d& N C NLO + ^NNLO^ = MnNLO [ d $m+l(P3, ■ ■ ■ ,Pm+3; X ± p l7 X 2 p 2 )-^- 
Jl J <->m- 



m+l Xi X 2 



1 

X - 

2 



^ik,j( s ik,xi,x 2 ) - X° kJ (s Jk -,x 1 ,X2) 



X° 3 ({p} m+3 ) 



X 



\M m+ 2({p}m+2)\ 2 jt\{p}m+2) , (2.73) 



together with correction terms coming from the oversubtraction of the single unresolved 
limits from the X® antenna (see Eq. 2.38) ultimately cancel against the integrated wide 
angle soft radiation term J x da^ e NL0 . 

To make this cancellation more precise, consider the coefficients of terms proportional 

to 

X? ek \M m+ 2(. ..,a,I,K,b,.. .)\ 2 Jt\{p}rn + 2)- (2.74) 

Here we consider final-final radiation, but the argument is general. The wide angle soft 
contribution is given by Eq. (2.43). As discussed in Section 2.1.3, the leading singularity 
from the wide angle soft terms is proportional to 

1 , ( SaiSl K Skb\ . 

— log +0(1). 

2e \SalSikSKbJ 

This pole cancels against combinations of the form (where the dependence on x±, X2 is 
suppressed), 

X - [X°(s tk ) - X°(s ai ) - Xl{s kb ) - Xl(s IK ) + Xl(s aI ) + X°(s Kb )] . (2.75) 

The first three terms in Eq. (2.75) are obtained by integrating terms in the double real 
subtraction term da^f^Q . The first term is produced by terms coming from f t d&^NLo 
that compensate for the oversubtraction of single unresolved poles associated with a four- 
parton antenna of the X® type as in Eq. (2.38). This can be understood as being due to 
repeated radiation where the hard radiators are particles i and k. The second and third 
terms have opposing signs (compared to the first term) and come from the almost colour- 
connected term j l da^ c NLO and occur when the integrated antenna describes unresolved 
radiation which took place between hard radiators a and i (b and k). For these three terms, 
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the unresolved radiation is emitted from an outer X® antenna, while X%, is the inner 
antenna. The final three terms form part of the real- virtual subtraction term da^jf LO , 
which we can understand as corresponding to emitting an unresolved momentum between 
pairs in the colour-connected set . . . ,a,I,K,b, . . . i.e. unresolved radiation from an inner 
antenna. In this case, X® ek plays the role of outer antenna. The signs are always fixed to 
be opposite of the partner term coming from J l d^vjvLO - As usual, the subtraction terms 
constituting d<7 nnlo mus t be integrated over the unresolved phase space and added back 
in integrated form to the double virtual contribution. Expanding (2.75), we find that the 
leading pole is proportional to 

1 , ( S a iSlKSkb\ . 
-7T log + 0(1 

2e \Sais ik s Kb J 

which cancels against the leading pole coming from the wide angle soft term. 
3. Renormalisation and mass factorisation 

We are concerned with m-jet production in the collision of two hadrons h±,h2 carrying 
momenta H\,H2 

M#l) + M#2) -> j(p?,) + ■ ■ .+j(Pm+2). (3.1) 

Within the framework of QCD factorisation the cross section for this process is written 
in Eq. (1.1) as an integral over the infrared and ultraviolet finite hard partonic scattering 
cross section for m-jet production from quarks and gluons, multiplied by parton distribu- 
tion functions (PDF's) describing the momentum distribution of these partons inside the 
colliding hadrons. 

However, after combining the real and virtual contributions together with the antenna 
subtraction terms, the partonic cross section dcx^- contains both ultraviolet and initial-state 
collinear singularities. We remove the UV singularities through coupling constant renor- 
malisation in the MS scheme and absorb the initial-state singularities into the PDFs using 
the MS factorisation scheme. For simplicity, we first set the factorisation and renormalisa- 
tion scales equal to a common scale, 

2 2 2 

V- = = Vf- 

It is straightforward to restore the dependence of the partonic cross section on these scales 
using the requirement that the hadronic cross section is independent of them. 

3.1 Ultraviolet renormalisation 

In terms of the bare (and dimensionful) coupling a b s , the unrenormalised cross section has 
the perturbative expansion, 

ckx™ = da^ LO + (£\ da™> NLO + (£\ 2 da™> NNLO + O ((a b s f) (3.2) 

bearing in mind that the LO cross section is O ((o^) m ). Furthermore, each power of the 
bare coupling is accompanied by a factor of (7(e) given by Eq. (2.11). 
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The renormalisation is carried out by replacing the bare coupling with the renor- 
malised coupling a s = a s {jjL 2 ) evaluated at the renormalisation scale ji 2 R = fj 2 , 



a° s = Z as (j%a s . 



In the MS scheme [92], 



1 



C(e) 

where /3q and f3\ are [93-96] 



1 



2tt 



+ 



e 2 2e 



a- \ 2 
2-kJ 



Y + 0(a 3 s ) 



(3.3) 



(3.4) 



/3i 



11C A 2Np 



UC\ 1QC A N F C F N F 



12 



12 



(3.5) 
(3.6) 



Substituting Eqs. (3.3) and (3.4) into Eq. (3.2) and collecting with respect to a S} gives 
the coefficients of the renormalised (but not infrared-finite) expansion, 



= *%> + (£) da?™ + (g) 2 d^ L ° + O (o») 



(3.7) 



.2irJ LJ \2ir 

Since the m-jet cross section has a leading order behaviour proportional to a™, we have, 



dor. 



LO 

ij 



da. 



un,LO 
ij 



da; 



NLO 

ij 



1 , un,NLO 

W) 11 



m/3o 



da™ 



un,LO 



da 



NNLO 



1 



un,NNLO 



(m + l)Po 1 un,NLO 

e W) 13 



mfii m(m + 1)(3 2 



2e 



2e 2 



i un,LO 



(3.8) 
(3.9) 

(3.10) 



Note that for the double real radiation contribution da§^ LO and its associated subtraction 
term da^ NLO , renormalisation simply amounts to the replacement C(e) a b s — > a s . Likewise, 
the inverse powers of C{e) are immediately cancelled against the additional factors of C{e) 

, ■ j un.NLO i , un,NNLO 

present in aa- and da- 
3.2 Mass factorisation 

After renormalisation, the physical cross section is given by 

#2 

• / ~ 

';3 



da(H 1: H 2 



6 



fii&fji&daij&H^&Hz) 



(3-11) 



where the bare PDF for a parton of type a carrying a momentum fraction £ is denoted by 
f a (0- For clarity, we have made the dependence of the partonic cross section ddij on the 
initial state parton momenta, £,\Hi and £2-^2) explicit. 
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The initial-state singularities present in dcr^ are removed using mass factorisation to 
produce the hadronic cross section where the finite partonic cross section multiplies the 
physical PDF / a (£, Mf) a * the factorisation scale \j? F = fi 2 R = /U 2 . The physical PDF are 
related to the bare PDF / a (£) by the convolution, 



(3.12) 



/ a (£, fJ 2 ) = J dxdy f b (x)T ba (y, /x 2 )<5(£ - xy) = [f b r 6a j (£, ^ 2 ) 
The kernel r ba has the (renormalised) perturbative expansion, 



T ba (y,n 2 ) = 5 ba 5{l-y) + 
where in the MS scheme, 

T ba(y) = --Pbaiv); 



2tt 



rL(y) + 



a s (p) 
2vr 



T 2 ba (y) + 0(a 3 s ) (3.13) 



r L(y) = ^2 



Z>k®*&J (y) + 2/3oA(y) 



(3.14) 
(3.15) 



where the are the standard four-dimensional LO and NLO Altarelli-Parisi kernels in 
the MS scheme given in Refs. [97-101] and are collected in the Appendix D for the gluonic 
channel. 

Eq. (3.12) can be systematically inverted such that the bare PDF is given by, 



f a (0 = I dxdy f b {x^ 2 )T^(y^ 2 )5(i-xy) = [h®T^] (£,^ 2 ) (3.16) 



where, 



r b - a W 2 ) = <WHi-y)-( 



2vr 



n a (y) 



2ir 

+0(al). 

Inserting Eq. (3.17) into (3.11) and applying the rescaling £j — > x^i we obtain 



tfa(y)-Y,l F lc® r ca\ (V) 

c 

(3.17) 



MHi,H 2 ) = W ^^fitiu^fji&^d&ij&H^bHz) (3.18) 

where the infrared-finite mass factorised partonic cross section is 

d& ij (Z 1 H 1 ,&H 2 ) = f ^ f ^Tll(x l ^ 2 )T7 l (x 2 ^ 2 )dd kl {x^ l H l ,x 2 ^H 2 ).{?>.l^) 
J Xl J x 2 J 

Expanding in the strong coupling as in Eq. (1.1), we find the connection between the 
infrared singular cross sections and the mass factorisation counterterms, 

d*ovLo(£itfi,6#2) = daij, LO (ti tfi,6# 2 ), (3.20) 
teij,NLo(SiH 1 ,Z2H 2 ) = dfiwLofaHi&Hz) + d^iofotfi, &H 2 ), (3.21) 
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da i j j NNLo(^iHi,^2H2) = da i j i NNLo(^iHi,^2H2) + da^Jf NLO (^iHi,^2H2), (3.22) 



where 



&Vij,NLo(£l H li& H 'z) = 



d&ij,NNLO (£1 H i ; 6 #2 ) 



' I -p Lr ki( x i) d Vkj,Lo(xitiHi,&H 2 ) 
— -rl j {x2)daii tL o(£iH 1 ,X2&H2), 

X2 

- [ —Tl i (x 1 )da kjtLO (xi(,iH 1 ,^2H 2 ) 

J X\ 

- I —T 1 ki (x 1 )da kjtNL o(x 1 ^ 1 H 1 ,^ 2 H2) 

J X\ 



(3.23) 



/ 
/ 



—tfj(,X2)dau,Lo(£iHi,X2&H2) 

—rl j (x2)da i i jNLO (£iH 1 ,X2&H2) 

^ I ^T 1 H {x l )T} J {x 2 )da^{x^ 1 H 1 ,X2^H2). 

X\ J X2 



Recalling that 
da ij:NLO = / (da^ 

Jd$ m+1 



NLO 



d&ij t NLO 



>+/ (/ 



d&ij t NLO 



(3.24) 



+ da(j NLO + dafjJ[ LO j , 
(3.25) 



we find that the NNLO mass factorisation counterterm for the (m + l)-parton phase space 
d °ij,NNLO in equation (1.4) is given by, 



d °ij,NNLO (tiHi ,&H 2 ) = 



-I 



dxi dx2 

X\ X2 

dx\ dx2 

X\ X2 



S(l-X2)T 1 ki (x 1 ] 

5(i- Xl )r} J (x 2 ) 



da 



R 



"kj,NLO ~ da%^ NLO 



da 



R 



l,NLO 



,NLO 



(x 1 £, 1 H 1 ,x 2 £,2H 2 ) 
{x 1 ( 1 H 1 ,x 2 ^2H 2 ), 



(3.26) 



while da^^ff LO in equation (1.5) contributes to the m-parton phase space counterterm, 



d ^ii,NNLO (£l H l, 6 H 2 ) 



+ 5(i - xi) % ( r* (x 2 ) - ^ [i% ® ry (* 2 ) 

-rL^i)^-^) f da kl>L o(xi€iH 1 ,X2bH 2 ) 
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5(1 - * 2 ) rL(^i) / da s kj:NLO + da 
j J-i J-2 L/i 



kj,NLO 



V 

iLNLO 



(x 1 £, 1 H 1 ,x 2 £,2H 2 ) 

{x 1 £ 1 H 1 ,x 2 &H 2 ), 

(3.27) 



and will be discussed elsewhere. 

3.3 Scale dependence of the partonic cross section 

We start from the hadronic cross section that depends on up = i~ir = [i through the strong 
coupling and the PDF's, 



h3 



The scale variation of the coupling constant a s (//) is given by, 

2 9 ( a s (/x) s 



(3.28) 



/x 



<9/i 2 



2vr 



= -/3(a s (/i)), 



where, 



£M/i)) = A) 



q s (m) 

2vr 



2vr 



(3.29) 



(3.30) 



with Po,Pi given in (3.5), (3.6). Similarly, the scale variation of the parton distribution 
function fi{x,y?) is determined by the DGLAP evolution equation, 







where, 



diJ? 



Pij(x) = p%{x) + 



2vr 



2vr 



[Pij ®fj{v 2 )] (x), 



(3.31) 



(3.32) 



Demanding the independence of the physical cross section in (3.28) on the unphysical 
scale fi 2 , 

.2 d 



:da = 0, 



(3.33) 



we obtain the following scale variation equation for the partonic cross section, 



(i 



_d_ 



2vr 



Pik ®da kj (a s (fi),fi) +da ik (a s (n),iJ,)®p kj 



(3.34) 



Bearing in mind that the leading order m-jet cross section depends on a™, and solving 
Eq. (3.34) order by order in a s , we see that the partonic cross section at scale fj, 2 is related 
to that at scale fii via, 
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daij(a 3 (n2),H2) = da ijiLO {a s (^ 2 ), Mi) 




ij,NLo{ot s {n2)^i) + m Pq Lda ij)L o{a s {ix 2 )^i) 



Lpki®da kjtL o(a s (fi 2 ),»i) - L p kj ® da iktLO (a s (n 2 ), Mi) 



+ 




^— J ^d<Ty )A rAr L o(a s (M2),Ml) 
+(m + 1) /3 £ da ijtN Lo(a s (fj,2), Mi) 



/ m(m + l) _o r2 \ „ / / \ \ 
+ lmpiL-| - p L j dai jjL o{a s (fi 2 ),fii) 




where L = ln^/Mi)- 

4. Real- Virtual corrections for gluon scattering at NNLO 

In this section we discuss the amplitudes that enter in the implementation of the mixed 
real-virtual correction. We focus on the pure gluon channel and describe the colour decom- 
position of the gluonic matrix elements at tree and loop-level. The remaining contributions 
that we use for subtraction, namely the three parton tree-level unintegrated and integrated 
antennae as well as one-loop three parton antennae have been derived in previous publica- 
tions. For convenience and completeness we collect them in appendix B. 

4.1 Gluonic amplitudes 

The leading colour contribution to the m- gluon ra-loop amplitude can be written as [78, 



where the permutation sum, Sjyi I 'Zjjxi, is tli© group of non-cyclic permutations of m symbols. 
We systematically extract a loop factor of (7(e)/2 per loop with (7(e) defined in (2.10). The 



102-105] 




a{m)). 



(4.1) 
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helicity information is not relevant to the discussion of the subtraction terms and from now 
on, we will systematically suppress the helicity labels. The T a are fundamental represen- 
tation SU(N) colour matrices, normalised such that Tr(T a T b ) = 5 ab /2. .A™(1, ■ ■ ■ ,n) 
denotes the n-loop colour ordered partial amplitude. It is gauge invariant, as well as being 
invariant under cyclic permutations of the gluons. For simplicity, we will frequently denote 
the momentum pj of gluon j by j. 

At leading colour, the tree-level (m + 2)-gluon contribution to the M-jet cross section 
is given by, 

da = A/^ +2 d$ m (p 3 , . . . ,Pm+2-,Pl,P2)^ 

x £ A° m+2 (a(l),...,a(m + 2))J^\p 3 ,...,p m+2 ) (4.2) 

where the sum runs over the group of non-cyclic permutations of n symbols and where 
A° m (a(l),...,a(m)) = £ A%(a(l), . . . , a(m))A° m (a(l), . . . , a(m)) + O (JL Y (4.3) 

helicities 

For tree-processes involving four- and five-gluons, there are no sub-leading colour contri- 
butions. The normalisation factor A/^ +2 includes the average over initial spins and colours 
and is given by, 

,ro ,r fa s N\ m ~ 2 C{e) m - 2 

^ + 2=^ox[-^) (4.4) 

with 

* r "> = h* 1)2 X (9 2 N) 2 (N 2 - 1), (4.5) 

and where we have absorbed the factors of g 2 using the useful factors C(e) (2.10) and C(e) 
(2.11), 

'a, AT 



9 NC{e) = \^r) c{e) - 

The leading colour six- gluon real-real contribution to the NNLO dijet cross section is 
obtained by setting m = 4 and M = 2 in (4.4) such that in this case J^^nlo appearing in 
Section 2 in (2.4) is given by N$§lo = ^6- 

For convenience, we introduce the additional notation for the one-loop "squared" ma- 
trix elements 

A 1 m (a(l),...,a(m)) = £ ^ A%(a(l), . . . ,a(m))A l m ^a(l), . . . ,a(m)) + O (^] 

helicities i=0,l 

(4.6) 

so that the one-loop, (m + 2)-gluon contribution to the M-jet cross section is given by, 
do- = A/^ +2 d$ m (p 3 , . . . ,p m +2-,Pi,P2)^ 
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x A 1 m+2 (a(l),...,a(m + 2))J^\p 3 ,...,p m+2 ). (4.7) 

As before, the normalisation factor includes the average over initial spins and colours 

and is given by, 



m+2 - MloX \^) cW^- (48) 



We will encounter both A\ and A\ when computing the real-virtual corrections relevant 
for the NNLO dijet cross section. 

4.2 The five-gluon real-virtual contribution da^ NLO 

The leading colour five-gluon real-virtual contribution to the NNLO dijet cross section is 
obtained by setting m = 3 and M = 2 in (4.8) such that in this case M^Yjlo = -^5 and 
we have, 



^NNLO =^L0d*3(P3,-..,P5;Pl,P2)^ £ A 5^(l), ■ ■ ■ M*))^ (P3, -.Ps) 

/a s iV\ 2 (5(e) 2 2 
= ° 1 ~27r~ J "^)" d 3 ^ 3 '---' P5;Pl ' P2 ^3! 

X XT I ^(Iff) iff, feff) + Al(lg, Ig, 2g, jg, kg) I Jg ^(Pi,Pj,Pk) 

P(i,j,k)e (3,4,5) V / 

(4.9) 

where the sums runs over the 3! permutations of the final state gluons. 

Therefore, depending on the position of the initial state gluons in the colour ordered 
matrix elements, we have two different topologies. These are labelled by the the colour 
ordering of initial and final state gluons. We denote the configurations where the two 
initial state gluons are colour-connected (i.e. adjacent) as IIFFF and those where the 
colour ordering allows one final state gluon to be sandwiched between the initial state 
gluons is denoted by IFIFF, 

A-RV A ~RV,IIFFF . A ~RV,IFIFF , ^ 

d °NNLO = d °NNLO + d<J NNLO > 



where, 



a ~rv,iifff Kr ( a s N \ 2 C{e) 2 

da NNLO =Nl0 \2^~) "^7gy d$ 3(P3,-..,P5;Pl,P2) 



X 3J Al(ig,2g,ig,jg,kg)4 3) (P3,---,P5), (4.11) 

P(jj,fe)e(3,4,5) 



A ~RV,IFIFF K r f a s N \ 2 C(e) 2 

da NNLO =Nl °\2^J -^- d ^{P3,---,P5-,Pl,P2) 
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E Al(i g ,i g ,2 g ,j g ,k g )4 3) (p 3 ,...,p 5 ), (4.12) 

' P(i,j,fc)G(3,4,5) 

The one- loop helicity amplitudes for gg — >• ggg have been available for some time [106]. 
We have cross checked our implementation of the one-loop helicity amplitudes of Ref. [106] 
against the numerical package NGluon [107] . 

We note that the renormalised singularity structure of the contribution in (4.9) can be 
easily written in terms of the tree- level squared matrix elements multiplied by combinations 
of the colour ordered infrared singularity operator [76] 

Therefore the real-virtual correction contains only jW type of operators and in the gluonic 
approximation I g l g is the only operator that appears. The singular part of the renormalised 
colour ordered gluonic amplitude takes the form, 

A-sO-g, ^gi^giJgi kg) = 

2 (l$(e,si 2 ) + I$(e, S2 i) + I$(e,*v) + I$(e,s jk ) + I$(e,s kl j) A° 5 (l g ,2 g ,i g ,j g ,k g ) + O(e ), 
A-50-g, *s> ^gi3g-> ^g) = 

2 (l$(e, s u ) + I«(e, s t2 ) + I«(e, s 2j ) + lg(e, Sj - fc ) + I$(e, Sfcl )) 2 9 , j 9 , fc fl ) + O(e ). 

(4.14) 

In Section 5 we will explicitly write down the counterterm that regularises the infrared 
divergences of the real-virtual correction for topology (4.11) and (4.12) separately. 



1 gg \ b gg) 



on 



2r(l-e) 



1 11 

+ 6e 



5. Construction of the NNLO real-virtual subtraction term 

As stated in the introduction, the aim of this paper is to construct the subtraction term for 
the real-virtual contribution such that the (m + l)-parton contribution to the m-jet rate 
is free from explicit e-poles over the whole of phase space and the subtracted integrand is 
well behaved in the single unresolved regions of phase space. As discussed earlier, this is 
achieved with the help of the antenna functions and, as will be explained here and in the 
subsequent section, the limit e — > can be safely taken and the finite remainders evaluated 
numerically in four dimensions. 

We start this section by recovering the general formula for the real virtual channel 
which has to be integrated over the (m + l)-parton final state phase space numerically. It 
reads, 

/ ( d&NNLO ~ d °NNLO + / d °NNLO + d °NNLO ) ' ( 5 - 1 ) 

Jd$ m+1 V Jl J 

where da^-^LO contains the once integrated part of the integrated double real subtraction 
term discussed in Section 2.1, da^ff LO is the the real-virtual subtraction term discussed 
in Section 2.2 and the mass factorisation term d&NNLO * s Si yen m Eq. (3.26). 
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The remaining contribution from the double real subtraction term denoted by J 2 da^ 2 NLO 
must be integrated over two unresolved particles and contributes directly to the m-jet fi- 
nal states. It must therefore be added to the integrated real-virtual subtraction term 
j 1 da^f%io ' t ne two- loop matrix elements ^Yinlo together with the mass- factorisation 
counterterm denoted by da NN ,2 Q . These all contribute to the m-parton final state and will 
be treated elsewhere. At the end of this section, however, we present the contribution from 
d&Yf%io wmcn must be added back in integrated form to the m-parton final state. 

As discussed in Section 2, the phase space for the (m + l)-parton phase space can be 
written as an integral over the longitudinal momentum fractions x±,x 2 in the form 



f d$ m+ i(p 3 , . . . ,p m +3;XiPi,X 2 p2)-^- 
J <->m4 



dx± dx2 



Wl x l x 2 

to account for initial state radiation. 

In the next subsections we present separately the subtraction terms required for gluon- 
gluon scattering divided for the IIFFF and IFIFF topologies. Subtraction terms related to 
each topology are denoted with an superscript X 5 and Y 5 respectively such that, 6 

d&NNLO = ^^NNLO + ^^NNLO' 0>-2) 

When presenting the subtraction term we group the terms which are free of explicit 
e-poles in square brackets. For conciseness we suppress the explicit x±,x 2 dependence of 
the integrated antennae appearing in the formulae below. 

5.1 IIFFF topology 

The one-loop single unresolved subtraction term for the IIFFF topology is: 
'a s N\ 2 C{e) 2 2 ^ f d Xl dx 2 
P(i,j,k)e(3,4,5)' 

l^l(s- 2i ) + \^( Sij ) + \j$(8 jk ) + ^W) + ^(«12) " ^(Xi)^! 



,. TZ , fa s N\ 2 C(e) 2 2 v-^ [ dx x dx 2 ^ , _ 

d °NNLO=*LO^) L J — — d*3(p3,..,P5;Pl,p- 2 )< 



X 2 ) 



-F 1 gg (x 2 )6(l - Xl)jA%(ig,2g,ig,jg,kg) J^fa , P J , P ' k) 

+fl(%, ig ,j g ) 5(1 - xi)<5(i - x 2 )A\(ig, %, (ij) g , k g ) + (rIKsn) + 

+\^ s Jif )k ) + \H{sik) - r 1 gg (x 1 )S(l - x 2 ) - T l gg {x 2 )5{l - x^ 



xA A (l g ,2 g ,(ij) g ,kg) 



J 2 \pij,Pk) 



6 This notation follows naturally from that used for the double real subtraction term [75] where the 
IIFFFF, IFIFFF and IFFIFF topologies were labelled Xe, Ye and Ze respectively, and anticipates a similar 
labelling in the double virtual contribution where the IIFF and IFIF topologies will be labelled X4 and 
Y4. We note that X5 receives contributions from singly integrated contributions from Xe and Ye, while Y5 
receives contributions from Ye and Z@. 
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+f3(igi3gi kg) 



5(1 - x,)5(l - x 2 )A\(t g , %, (ij) g , (jk) g ) + ( J^(s V2 ) + ^!(s m ) 



3" 6 y (ij)Uk) 
xA 4 (lg,%,(i3)g,Gk)g 



+ - V gg^) 5 ^ - x *) - T lg( x 2)S(l - x^ 

J?\pij,Pjk) 



+f3{jg,k g ,lg) 



5(1 - Xl )S(l - X 2 )A\(lg, %,i g , (kj)g) + ( J^(s^) + -Ii(s 2i ) 



XA^lg, 2g,ig, (kj)g) 



J 2 \pi,Pkj) 



+ 



fl(%, i g ,j g )6(l - Xl )5(l - x 2 ) + ( \j^(s 2i ) + l -H(s 2j ) + l -Fl( Sl] ) - J^(s m ) 



X fi(2g,ig,jg) 



Al(lg,2g, (ij)g, kg) (pij,P k ) 



+ 



fl(tg,jg, kg)5(l - Xl )5(l - X 2 ) + ^ (S tj ) + ^(S ik ) + ^(Sj,) " ^(^^ 



x f$ (igijgi kg) 



Al(l g , 2 g , (ij)g, (jk) g ) 4 2) (pij,p jk ) 



+ 



fl(j g , kg, l g )8(l - Xl )5(l -x 2 ) + [ ^(s jk ) + \ft( Slj ) + l -Fl(s lk ) - Fl(s 



x /3 {jgi kg, lg) 



Ai(lg,2g,ig,(kj)g)4 2) (pi, Pkj ) 



11 , 

+ ylog 

11 . 
+ylog 

11 . 

+7r lo s 



1 

+ 2 



f* \ to 



>2ij\ 



fi(%, ig,jg)6(l ~ Xl )5(l ~ X 2 ) A 4 (lg, = 2g, (ij)g, kg) fa , Pk ) 



fsiigJg'kg) 6 ^ ~ ^i)^ 1 ~ x ?) A l0-g> 2 g> U k )g) J 2 \PijiPjk) 



» \ to 



hjk\ 



/s k g> 1 g) 6 ( 1 ~ xi)5(l - x 2 ) A° 4 (l g , 2g,i g , (kj)g) jf ] (pi,p kj ) 



1 



\^{s- 2W) ) + l -Fl(s- 2l ) - \n{s m ) + \^(s lk ) + \n{s mm ) ~ ^(-*) 



2(ij),ik) 5 ( S 2j ) s ik,X 2 i,i k ) + S(s^. k yS ik , X l( j k ^ ik ) 



+5(1 - Xl )5(l - x 2 ) (S(s^,s ik , 
S(si k , s ik , xi kyik ) - 5(s^ ( ~ } , s ik ,x^^ ) . k ) + S(s ik , s ik , 1)) 

X f3(igJgikg)A 4 > (tg,2 1 g, (ij)g, (jk)g) ^{PijlPjk) 
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1 

+ 2 



^(Slj)- ^3(^12) +^( S 12) 

+<Hl-x 1 )<5(l-x 2 )(s( S . ( ~ ) ,^^ 



/3 (Ifl) kg,jg) Al(lg, 2g, lg, (kj)g) \pi,ftkj) 



1 

+ 2 



1 



+5(1 - ZiWl - X 2 ) ( S{s k ^ Y S kj ,X kJ rr )k .) - S(s kj , S kj , 1) - S{s=jrr v S kj , X=^ fc . 



k(ij),kj 



2(ijy"K3'"2(ij),kj> 



+S(s 2j ,S kj ,X 2jykj )^j 



/3 U ( 2 9 , *» i ) ^4 ( Iff . 2 9 . (U )<? , kg ) J 2 (2) (ftj , Pk ) 



1 

+ 2 



^"3(^12) + ^12) - ^3 (%■) + \H(s2j) - \^ik) + \^l(nk) 

+5(1 - xi)5{l - x 2 ) (S(sis, Sj- k , x i5 y k ) - <S(s l2 , Sj k , + <S(%-, s jfc , x 2jJjfc ) 

-<S(s§j, s^, x^.jji) + 5 ( s ifc> s jk,xi kjjk ) - 5 ( s Ifc> s jfc' 

X (l fl , ig , 2 9 ) Al ( 1 9 , 2 9 , jg , kg ) , ^ ) ^ . 



5.2 IFIFF topology 

The one-loop single unresolved subtraction term for the IFIFF topology is: 



d<r T ' Y5 

NNLO 



fa s N\ 2 C(e) 2 2 ^ fd Xl dx 2 ^. _ _ 

Mlo \-*T) W)^. ^ y— — d*3( P3 ,...,P5;Pi,P 2 )« 

V 7 W P(iJ,fc)6(3,4,5) 

^(-li) + + + ^(-i*) + l^kl) - Tig(xi)S{l - x 2 ) 

-r gg (x 2 )<5(l - xi) )Ag(i ff ,i ff ,5 9 ,j ff ,fc ff )^ 3) (pi,pj,p fc ) 



+F 3 °(l 9 ,i 9 ,2 9 ) 



5(1 - sxWl - x 2 )A|(l s ,2 ff , j g ,~kg) + (^(512) + -J^) 



+^(«5*) + ^§(*is) - ^(^i)^ 1 - **) - r 99 (x 2 )5(i - xx)) 



x^-4(lff) ^giJgi kg) 



J?\pj,Pk) 



+fii^gi3gi kg) 



5(1 - Xl )5(l - X 2 )A\(ig,ig,%, (jk)g) + ( ^(Sfc) + ^(%) 



(5.3) 
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'"A^lgjlg, 2g, (jk)g) 



J 2 \pi,Pjk) 



S( I ,r i)( )( ! - .r 2 )A\a ir /„. %,. (}7j) fl ) + ( ^(.s u ) + \^{s 2l ) 



+^( s 2(1yj) + ^(*i ( Ty)) " ^(xi^l - x 2 ) - r^(x 2 )<5(l - x0) 

J^HPi'Pkj) 



+ 



x^-4(l g ,i g ,2 g , (kj)g) 



^ ( ','/• V 2,, W 1 -n I .r 2 ) - ( Ijf + Fi(s r2 ) + lj? ( S2i ) - 2F$(s 



12- 



xF 3 (l 9 ,i g ,2 g ) 



A2(l ff ,2 fl ,j ff ,fe p ) 4 2) {pj,Pk) 



+ 



fffigJg, k g )6(l - Xl )6(l - x 2 ) + [ \^{s 2j ) + \j^{s 2k ) + ^(s jk ) - J*(s 



x fi(^gi3ai kg) 



2 ov ZJ/ 2 



+ 



1 



1 



fUjg, kg, i g )5(l - Xl )6(l -x 2 )+[ ^(s jk ) + -H{ Slj ) + -Tl(s- lk ) - J^(s m ) 



x fi (jgi kg, Iff) 



A2(l fl , i fl , 2 fl , J 2 (2) {Pi,p k j) 



11 



+ Y ( |^T|) F i(%ih\)^ -Xl)5{l -X 2 )Al(lg\,jg,~kg) jf , Pfc) 



11 



+ y lOg [j^) f!(%dg,kg)5(l - Xl )6(l - X 2 )Al(lg,ig,2g, (jk)g) J^ijU^jk) 

+11 i g f t Q (A u 1 wi _ ^.uh _ „„\ ^o^T „■ o o^h ^ 7-( 2 )/ 

6 Vl*ijfcl. 



r(2)< 



fi{jg,kg,l g )8(l -a:i)<5(l - ^2)^4(15, i 9 , 2 9 , 4 (PhPkj) 



+ 



+- 



7j(a n ) - J^( Sn ) + \j^{s 2j ) - lj%{aq) + ^fe) - ^(s Vk ) 
+5(1 - xi)8(l - x 2 ) ( - 5(sj2, sjfc, a?Hjfc) + 5(s l2 , s jjfc , z l2 - Jjfc ) - 5(s 2j , s jfc , x 2jJk ) 
+S(s^,s- jk ,x rj - jk ) - S(s lk ,s jk ,xi kJk ) +S(si k ,s- jk ,xi k - jk ) S j 

xFl{\ g ,igX)Al{%X^k)J?\p^Pk) 
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s 2(kjy Si ^ x 2(kj),ij 



) ^( s 2jr'> Sij,X 2 j t ij) + Sjj, Xijjj) 



)) /3 (Iff' is) ^(Ig^gi 2- 



(kj) g )4 '{PhPkj) 




1 



" k 3 ° (s iw ) + ^ (aifc) " ^ (s2fc) + k 3 ° ( %5)) " ^°(^) + ^(%) 



+5(1 - - X 2 )^(s I( ~ ) ,S ife ,X I( ~ ) . jfc ) -S(si k ,S ik ,Xl k}ik ) +S(S2 k ,S ik ,X2k !ik ) 

~ ^( S 5(fci)' S * fc ' X 2(fcjj,ifc)) ^3 (^S'is' ^s) ^4(^S^S> ^35 (i^) ff ) ^2 \PiiPjk) \- 



(5.4) 



For this topology, the contribution to the hard region (xi,X2 7^ 1) turns out to be 
identically zero. This is because the only terms that contribute in the hard region are inte- 
grated initial-initial antennae. This configuration (4.12) does not involve colour-connected 
initial state gluons and therefore does not contain any integrated initial-initial antennae. 7 

5.3 Infrared structure 

With the explicit expressions for the integrated antenna functions given in Appendix B, 
the integrated large angle soft terms given in Section 2.1.3 and the pole structure of the 
one-loop matrix elements, it is straightforward to check that the explicit e poles analytically 
cancel in each and every one of the groups of terms in square brackets in Eqs. (5.3)-(5.4). 
The only remaining poles lie in the second lines of Eqs. (5.3) and (5.4) which precisely 
reproduce those of the real-virtual matrix elements. The limit e — > can therefore be 
safely taken. 

In summary, we have shown that, within the antenna subtraction formalism, the real- 
virtual corrections to gluon-gluon scattering are locally free of explicit e-poles providing us 
with a stringent check on the construction of the necessary subtraction terms. This is in 
direct contradiction to the statement made in Rcf. [32]. 

Furthermore, Eqs. (5.3) and (5.4) are free of implicit kinematical singularities corre- 
sponding to the single unresolved regions of the phase space as will be demonstrated in 
Section 6. 

5.4 Contributions to the m-jet final state 

In this subsection we identify the contributions from the real-virtual channel that we have 
subtracted in unintegrated form and which therefore must be added back in integrated 
form in the double virtual (m + 2)-parton channel. As expected contributions from the A5 
topology collapse in integrated form to the A4 topology of the virtual-virtual contribution 
while contributions from the Y5 topology contribute to both X4 and Y4. The contributions 

7 Due to the grouping of terms in Eqs. (5.4) to demonstrate the cancellation of e-poles, this type of 
initial-initial antenna does appears. However, they always cancel pairwise amongst themselves. 
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of da)^fj L Q, which when integrated over the antenna phase space become proportional to 
the X4 topology are denoted by dd'YfNLoIx an d are given by, 



da 



vs I 

NNLO\Xl 



K . (a s N\ 2 C(e) 2 2 ^ 

°l^rj ~Cie)^\ ^ d<$> 3 {p 3 ,...,p 5 ;p l ,p 2 ) 



P(i,j,fc)e(3,4,5) 



F 3 °(i ff , i fl , 2 fl )4(i fl , 2 fl , J 2 (2) (pj,p fc ) 

+/3 (2 9 ,i 9 ,i 9 )Al(i 9 ,2 9 ,(ii) 9 ,/c 9 ) 4 2) (pij,p k ) 

+/s(*ff> i»> M^K^ 2 9 , (tj) fl , 4 2) {pij,p jk ) 
+f$tig, k g , W A \6-g, 2 g ,i h , (kj) g ) J^iPuPkj) 

+ ^(Ifl.ifl,^) + ^F 3 °(l 9 , V 2,) (( Sl2j )- £ - (/x 2 )- 6 ) ) A°(i 9 ,2 9 , fe,p fc ) 



11 



+ ( ti(2 g ,i g ,j g ) + -/ 3 °(2,, V j 9 ) ((| S2 y|)- £ - (M 2 )" £ ) ) 4>(i fl ,2 fl , 4 2) (Pij,p k ) 



11 



6e 
11 



+ ( fhigJM + £f$(i g ,jg,k g ) ((s ijk )-* - (/x 2 )- 6 ) ) A°(i g ,2 9 , (ij) (jk) g ) jf\p l3 ,p ]k ) 



+ ( /sO,. l g ) + -/ 3 °(j 9 , fc ff , l g ) {{\s lkj \r ~ (V 2 r e ) ) A%(l g , 2 g ,i g , (kj) g ) 4 2) (pi,p kj ) 



r( 2 )i 



fa 3 N\ 2 C(e) 2 2 ^ /" da* da* , _ _* 



+F°(l g ,i g ,2 g ) 



P(i,j,k)e(3,4,5) 



: Al(lg,2g,jg,kg)4 ^(PjlPk) 



lk) 



+f3@g,ig,3g) 



:^4(1 9 > 2 9> (*j')«>, M 4 2) (.Pij,Pk) 



+f$ ^gijgi kg) 



x^(i 9 ,2 9 , (y) (jk) j 4 2) (Pij,Pjk) 



+f!i(jg,k g ,lg) 



xA2(l 9 , 2 ff ,i s , (/cj) s ) 4 2) (Pi'Pkj)- 
Similarly, the contributions which become proportional to the I4 topology are, 

ir /a s iV\ 2 (5(e) 2 2 ^ .f 
dcJ 7VJVLoly 4 ° =A/lo ^"^J CXe)"3! ^ d$ s(P3, • • • ,P5;Pi,P2)< 



(5.5) 



P(jj,fe)e(3,4,5) 



/3 U (2 P , jg , kg ) A\ ( 1 fl , i fl , 2 g , (j fc) ) ) (pi , p (j fc) ) 
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11 



6e 
11 



+ (/3 1 (2 g ,i 9 ,A ;g ) + ^/3 (2 g ,i g ,fc g )((| S2jfe |)- e -(/x 2 )- e )) A0(i g ,z g ,2 fl ,(ifc) g )J 2 (2) ( Pi ,p . fc) ) 



+ ( /s 1 ^, + ^/ 3 °&. i 9 ) {(\sikj\r - c^n ) ^(i g ,ig, 2 g , (kj) g ) 4 2) ( Pi ,p m ) 



r( 2 )/ 



2-7T / C(e) 3! , / xi X2 

7 V ' P(»J,fc)€(3,4,5) 17 



d$3(P3, • • • ,P5',Pl,P2)< 



2" o \- -11/ ' ^- O ^ 

x^(l g ,i g ,2 g , (j/c) ) J%(pi,p jk ) 



+f?,{jg, k gi l g) 



XA2(lg,Zg,2 g , (/Cj) ) jl{Pi,Pkj)- 



(5.6) 



6. Numerical results 

In this section we will test how well the real virtual subtraction term daJf NLO derived in the 
previous section approaches the real virtual contribution da^ NLO m au single unresolved 
regions of the phase space so that their difference can be integrated numerically in four 
dimensions. We will do this by generating a series of phase space points using RAMBO [108] 
that approach a given single unresolved limit. For each generated point we compute the 
ratio of the finite parts of da^ NLO and daJf NLO , 



R 



Finite(da^ v NLO ) 
Tinite(da% NLO ) 



(6.1) 



Here dcrjy^o ^ s the interference between the one-loop and tree-level five gluon matrix 
elements given by Eq. (4.9), and daJj NLO is the real virtual subtraction term given by 
Eqs. (5.3) and (5.4). The ratio R should approach unity as we get closer to any singularity 
showing that the subtraction captures the infrared singularity structure of the real virtual 
contribution. 

For each unresolved configuration, we will define a variable that controls how we ap- 
proach the singularity subject to the requirement that there are at least two jets in the 
final state with pr > 50 GeV where the jets have been clustered with the anti-A^ jet algo- 
rithm [109, 110] with radius R=0.4. The partonic center-of-mass energy y/s is fixed to be 
1000 GeV. 



6.1 Soft limit 

To probe the soft regions of the phase space, we generate an event configuration with a 
soft final state gluon k by making invariant Sij close to the full center of mass energy 
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Figure 2: (a) Example configuration of a single soft event with Sj,- 
R for 10000 single soft phase space points. 



s i2 = s - (b) Distribution of 



s\2- This kinematic configuration is pictured in Fig. 2(a). We define the small parameter 
x = (s — Sij)/s and show the distributions of the ratio between the real- virtual matrix 
element and the subtraction term for x = 1CP 5 (green), x = 10~ 6 (blue) and x = 1CP 7 
(red) in Fig. 2(b) using 10000 phase space points. The plot also shows the number of 
points that lie outside the range of the histogram. We see that the subtraction term 
rapidly converges to the matrix element as we approach the single soft limit. 

6.2 Collinear limit 
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Figure 3: (a) Example configuration of a single collinear event with Sij 
R for 10000 single collinear phase space points. 



0. (b) Distribution of 



Next we probe the final and initial state single collinear regions of the phase space. 
These event topologies are depicted in Figs. 3(a) when gluons i and j become collinear, 
and 4(a) where gluon i becomes collinear with the incoming gluon 1. 

For the final-final collinear singularity, we introduce the small parameter x = Sij/s\2- 
Fig. 3(b) shows the distribution in R obtained for 10000 phase space points for x = 10~ 7 
(green), x = 10~ 8 (blue) and x = 10~ 9 (red). 
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Figure 4: (a) Example configuration of a single collinear event with su — > 0. (b) Distribution of 
R for 10000 single collinear phase space points. 

Similarly in the initial-final collinear limit, the small parameter is x = su/s\2 an d 
Fig. 4(b) shows the distributions of R for the same x- values of x = 10~ 7 (green), x = 10 -8 
(blue) and x = 10 -9 (red). 

As the small parameter x gets smaller, we see a systematic improvement in the conver- 
gence of the real-virtual matrix elements and the subtraction term. This is in contrast with 
the collinear limit of the double real all gluon subtraction term [75,111], but not surprising 
due to the simplicity of the final state where the partons are fixed to be in back-to-back 
pairs as shown in Figs. 3(a) and 4(a). Nevertheless, Figs. 3(b) and 4(b) show that the 
subtraction term does not approximate the real-virtual matrix element as well as in the 
soft limit (Fig. 2(b)). This is due to the presence of angular correlations in the matrix 
elements stemming from gluon splittings g —> gg. The collinear limits of tree and one-loop 
matrix elements are controlled by the unpolarised Altarelli-Parisi splitting functions which 
explicitly depend on the transverse momentum k± of the collinear gluons with respect to 
the collinear direction and on the helicity of the parent parton. As a result of this, the 
splitting functions produce spin correlations with respect to the directions of other mo- 
menta in the matrix element besides the momenta becoming collinear. These azimuthal 
terms coming from the single collinear limits vanish after integration over the azimuthal 
angle of the collinear system. This occurs globally after an azimuthal integration over the 
unresolved phase space. Here we are performing a point-by-point analysis on the integrand 
defined by the real-virtual matrix element and the subtraction term and because we use 
spin-averaged antenna functions to subtract the collinear singularities, the azimuthal an- 
gular terms produced by the spin correlations are simply not accounted for in the antenna 
subtraction procedure. 

To improve on this, several approaches have been discussed in the past. One possible 
strategy discussed in [21] is to proceed with a tensorial reconstruction of the angular terms 
within the antenna subtraction terms. A second approach is to cancel the angular terms 
by combining phase space points which are related by rotating the collinear partons by an 
angle of ir/2 around the collinear parton direction [58,112]. In this case, the azimuthal cor- 
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Figure 5: Distribution of R for 10000 single collinear phase space point pairs where the pair of 
phase space points is related by an azimuthal rotation of tt/2 about the collinear direction, (a) 
Final state collinear singularity (b) Initial state collinear singularity 



relations present in the matrix element at the rotated point cancel precisely the azimuthal 
correlations of the un-rotated point. 

This second procedure was demonstrated to be extremely powerful in improving the 
convergence of the double-real radiation subtraction contribution to dijet production in 
[75,111] in the pure gluonic channel. The strategy of combining pairs of phase space points 
related by a tt/2 rotation eliminated the correlations from both: the real-radiation and its 
subtraction term. In the latter case, the four-parton antennae are responsible for angular 
correlations. In the real-virtual contribution discussed in this paper, the correlations can 
arise in the real-virtual matrix elements &&nnlo an< ^ m ^ ne tree-level five gluon matrix 
elements present in the subtraction term. There is no contribution from the three-parton 
antennae. Therefore, the azimuthal effect is expected to be smaller than in the double-real 
case [75]. Looking at the distributions shown in Fig. 3(b) for x = 10~ 7 (green), x = 10 -8 
(blue) and x = 1CP 9 (red) and Fig. 4(b) for x = 1CP 6 (green), x = 10 -7 (blue) and x = 10 -8 
(red), we see that the correlations are clearly visible, but are indeed relatively small. 

Nevertheless, to eliminate the remaining azimuthal correlations, we show the effect 
of combining related phase space points discussed above in Fig. 5(a) for the final-state 
collinear singularity and Fig. 5(b) for the initial-state collinear singularity for the same 
values of the small parameter as in Figs. 3(b) and 4(b) respectively. We observe a sig- 
nificant improvement in the convergence of the subtraction term, particularly in the case 
of the initial-final collinear limit. The conclusion is that by combining azimuthally corre- 
lated phase space points, the antenna subtraction term correctly subtracts the azimuthally 
enhanced terms in a point-by-point manner. 



7. Conclusions 

In this paper, we have generalised the antenna subtraction method for the calculation 
of higher order QCD corrections to derive the real-virtual subtraction term for exclusive 
collider observables for situations with partons in the initial state to NNLO. We focussed 
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particular attention on the application of the antenna subtraction formalism to construct 
the subtraction term relevant for the gluonic real-virtual contribution to dijet production. 
The gluon scattering channel is expected to be the dominant contribution at NNLO. The 
subtraction term includes a mixture of integrated and unintegrated tree- and one-loop 
three-parton antennae functions in final-final, initial-final and initial-initial configurations. 
We note that the subtraction terms for processes involving quarks, as required for dijet or 
vector boson plus jet processes, will make use of the same types of antenna building blocks 
as those discussed here. 

By construction the counterterm removes the explicit infrared poles present on the 
one-loop amplitude, as well as the implicit singularities that occur in the soft and collinear 
limits. The e-poles present in the real-virtual contribution are analytically cancelled by the 
e-poles in the subtraction term rendering the real-virtual contribution locally finite over 
the whole of phase space. 

We tested that our numerical implementation of the antenna subtraction term be- 
haves in the expected way by comparing the behaviour of the finite parts of the one-loop 
real-virtual contribution dcj^^ rLO with the finite part of the real-virtual subtraction term 
daJf NLO for the five-gluon process in the regions of phase space where one particle is 
unresolved. The numerical convergence of these terms gives confidence that the infrared 
structure of the real- virtual matrix element is captured by antenna subtraction method in 
a systematic and accurate manner. 

The real- virtual subtraction terms presented here provide a major step towards the 
NNLO evaluation of the dijet observables at hadron colliders. Future steps include; 

(i) completion of the analytic integration of the initial-initial four-parton antennae. 

(ii) analytic cancellation of infrared poles between the analytically integrated antennae 
present in the subtraction terms and the two- loop four-gluon matrix elements. 

(iii) full parton- level Monte Carlo implementation of the finite four-, five- and six-gluon 
channels. 

(iv) the construction of similar subtraction terms, etc., for processes involving quarks. 

The final goal is the construction of a numerical program to compute the NNLO QCD 
corrections to dijet production in hadron-hadron collisions. 

Acknowledgment s 

We thank James Currie and Thomas Gehrmann for stimulating discussions. We thank Si- 
mon Badger for assistance in comparing with the NGluon package and Stephan Buehler and 
Claude Duhr for help in using the CHAPLIN package. This research was supported in part 
by the UK Science and Technology Facilities Council, in part by the Swiss National Science 
Foundation (SNF) under contract PP0022-118864, in part by an ETH-Fellowship program 
under grant No. FEL-15 10-2, in part by the Research Executive Agency (REA) of the Eu- 
ropean Union under the Grant Agreement number PITN-GA-2010-264564 (LHCPhenoNet) 



-46- 



and in part by the US National Science Foundation under Grant No. NSF PHY05-51164. 
EWNG gratefully acknowledges the support of the Wolfson Foundation and the Royal So- 
ciety and thanks the Institute for Theoretical Physics at the ETH for its kind hospitality 
during the completion of this work. AG and EWNG also thank the Kavi Institute for 
Theoretical Physics, Santa Barbara for its kind hospitality during the early stages of this 
work. 

A. Momentum mappings 

The NNLO corrections to an m-jet final state receive contributions from processes with 
different numbers of final state particles. In the antenna subtraction scheme, one is re- 
placing antennae consisting of two hard radiators plus unresolved particles with two new 
hard radiators. A key element of the antenna subtraction scheme is the factorisation of 
the matrix elements and phase space in the singular limits where one or more particles are 
unresolved. This factorisation is guaranteed by the momentum mapping. 

In this section we denote the set of momenta for the M-particle process by {p}m- 
In order to subtract a particular singular configuration in a given process, we derive sub- 
traction terms which reproduce the exact singular behaviour of the matrix element in the 
unresolved configuration and employ a momentum mapping to implement momentum con- 
servation away from the unresolved limit. This has the consequence of mapping a singular 
configuration in an M-particle process to an (M — 1) or (M — 2)-particle process, depending 
whether the given singular configuration involves a single or a double unresolved limit. In 
integrated form these subtraction terms have explicit e-poles and contribute to final sates 
with fewer particles. The consistent momentum maps we require are, 



Let us consider the single unresolved emission that is relevant in this paper - either 
as part of the integrated single unresolved subtraction term that cascades down from the 
double real emission process (A.l) or in the single unresolved limit of the real- virtual 
contribution (A. 3). If the antenna consists of an unresolved particle j colour linked to two 
hard radiators i and k, then the mapping must produce two new hard radiators / and 
K. Each mapping must conserve four- momentum and maintain the on-shellness of the 
particles involved. There are three distinct cases, 



where, as usual, initial state particles are denoted by a hat. In principle, the momenta 
not involved in the antenna are also affected by the mapping. For the final-final and 



Mm+4 -> Mm+3, 
Mm+4 -> Mm+2, 
Mm+3 -> Mm+2- 



(A.l) 
(A.2) 
(A.3) 



final — final configuration 
initial — final configuration 
initial — initial configuration 



ijk ->■ IK 
ijk ->■ IK 
ijk ->■ IK 
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initial-final maps, this is trivial. Only in the initial-initial case are the spectator momenta 
actually modified. The momentum transformations for these three mappings are described 
in refs. [71, 113, 114] and will be recalled below. 

A.l Final-Final mapping 

The final-final mapping is given in [113] and reads 

it U U i U i u 

Pi = v—^ = xPi + r Pj + z P k 

Pk = P~ = (l-*K + (l-r)p£ + (l-*)p£ (A.4) 



where, 



(1 + p) s ijk -2rs jk 



z = 



2(Sij + Sjfc) - 

9( . \, \ [C 1 - p) s w ~ 2 r s a 

2{Sj k + s ik ) L 



p 2 = l + 4r(1 - r)g ^ fc . (A.5) 

The parameter r can be chosen conveniently [113, 114] and we use r = Sj k /(sij + Sj k ). 

The mapping (A.4) implements momentum conservation prrr, + prrr, = Pi + Pj + p k 
and satisfies the following properties: 

"k=°- "k =0 - 

Pjjjj -> P ( 7^ ->• Pfc when i is soft, 

p^rr — >■ pi + pj , VTu\ ~~ ^ Pk when i becomes collinear with j , 

^W) ~~ ^ ~* ^ ^ k w ^ en j becomes collinear with k. 

A. 2 Initial-Final mapping 

The initial-final mapping is given in [71] and reads 

It -U - it 

Pj = Pi = XiPi, 

P» K = = $ + Pi ~ (1 - *0 Pf , (A.6) 
with = p 2 K = and where Xi is given by, 

Xi = s Ji±Hh±l2lL . (A.7) 

Sjji ~\~ Si k 

Proper subtraction of infrared singularities requires that the momenta mapping satisfies, 

XiPi Pi , PTu\ ~~ ^ Pk when j becomes soft , 

XiPi Pi , P(jk) ~~ ^ P? ^ P k when j becomes collinear with k , 

XiPi p%— Pj , P(jk) ~~ ^ P k when j becomes collinear with i . 
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A. 3 Initial-Initial mapping 

The initial-initial mapping for ik,j — > IK is given in Ref. [71] and reads 







= XiPi, 












= Pe - 



{q + q) 2 q 2 
where p 2 . = p 2 ^ = p 2 = and 

T = Pi+Pk- Pj, T = Pi+Pk ■ 

The rescaling of the initial state momenta are given by the fractions Xi and Xk which 
read [71], 



Xi 



Sik + Sjk 


j Sik + Sij + Sjk 


Sik + &ij V 


' s ik 


Sik + Sij 


j Sik + Sjk + Sij 


Sik + Sjk \ 


' s ik 



Xk = . r 1 *^ 3 ^^ , (A.9) 

V Sik + s jk y s ik 

With these definitions it is straightforward to check that the momenta mapping satisfies 
the following limits required for proper subtraction of infrared singularities, 

Pi Pi , Pk — ^ Pk when j becomes soft , 
pi — >■ (1 — Zi)pi , pk Pk when j becomes collinear with i , 
pk — > (1 — Zk)Pk , Pi — ^ Pi when j becomes collinear with k . 

B. Gluonic antennae 

In this section we collect the three-parton tree-level unintegrated and integrated antennae 
as well as one-loop three-parton antennae used in the implementation of the real virtual 
subtraction term of Section 5. 

B.l Tree- level three-parton antennae 

The tree-level antenna functions are obtained by normalising the colour-ordered three- 
parton tree-level squared matrix elements to the squared matrix element for the basic 
two-parton process, 

yO _ q \-Mjjk\ 
^ijk — D ijk,IK ^0 |2 ' 

where S denotes the symmetry factor associated with the antenna, which accounts both for 
potential identical particle symmetries and for the presence of more than one antenna in the 
basic two-parton process. For the purposes of this paper the relevant antennae are, those 
containing the pure gluon final state. They are obtained from the radiative corrections to 
Higgs boson decay into gluons [57]. 
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B.l.l Final-Final emitters 

The tree-level three-parton antenna corresponding to the gluon-gluon-gluon final state 
is [21]: 

F 3 °(l 2 3 ) = — | g l 23g12 + S 123 g 13 + ^123 fi 23 + ^12513 + ^12^23 + S n S 23 
S 123 \ S 13S23 S 12 S 2 3 Sl2«13 S 23 «13 «12 

+4s 123 + 0(e) , (B.l) 

where Sij = (pi + Pj) 2 - As can be seen from the pole structure, this tree-level antenna 
function contains three antenna configurations, corresponding to the three possible config- 
urations of emitting a gluon in between a gluon pair. We make the decomposition [21] 

F 3 °(l, 2, 3) = / 3 °(1, 3, 2) + / 3 °(3, 2, 1) + / 3 °(2, 1, 3), (B.2) 

where 

/ 3 °(1, 3, 2) = f 2 fk£i£ + flHfH + £12£23 + 8 + \ 

S 123 \ • S 13 s 23 -S23 S 13 3 / 

The subantenna / 3 (i, j, k) has the full j soft limit and part of the i \\ j and j \\ k limits of the 
full antenna (B.l), such that i and k can be identified as hard radiators. Each subantenna 
is associated with a unique {3— >2} momenta mapping, k) — > (I,K) [113,114] given by 
Eq. (A.4). 

B.1.2 Initial-Final emitters 

The initial-final three-gluon antenna function can be obtained from its final-final counter- 
part (B.l) by the appropriate crossing of one of the particles from the final to the initial 
state, i.e. by making the replacements, s 2 3 ->■ (P2 + P3) 2 , S12 -> (Pi -P2) 2 , «13 -> (Pi -P3) 2 
and Q 2 = s 12 + S13 + S23- It reads [71], 

pO,i 9 o ^ 1 ( H2 , H2 , 8s? 3 , 8^?3 , 8s| 3 , 8s| 3 

^3 l A 9) l 9-> 6 9) — 0/n2 \2 1 1 1 1 1 

2(Q 2 Y y S 13 S 23 S12 S 23 S12 S13 

12si2Si3 12g 23 g 13 I2S12S23 4sf 2 4sf 3 4s| 3 



S23 «12 «13 S23S13 S23S12 S12S13 

+24s 23 + 24s 12 + 24s 13 + O(e) , (B.4) 

where the hat identifies the gluon crossed to the initial state. Because the initial gluon can 
never be soft, it is convenient to decompose this antenna into two contributions, that each 
contain the soft singularities of one of the final state gluons, 

F 3 °(l, 2, 3) = / 3 °(i, 2, 3) + / 3 °(1, 3, 2). (B.5) 
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Here, 



/3°(U3) = + ^ + + 4S - 



2(Q 2 ) 2 V S 12 S12 S12 S23S12 S 12 (S12 + S 13 ) 



| 8sf 3 , 6si 2 si3 



+ 512513 + 12s 23 + 12si2 + 12si 3 + 0(e) ) . (B.6) 

«23 S23 ' 

The subantenna f®(l,j, k) has the full j soft limit, the full 1 || j limit and part of the j \\ k 
limit of the full antenna (B.4), such that we can identify 1 as the initial state radiator 
and k the final state radiator. To numerically implement this antenna we use the {3—7-2} 
mapping, (i,j,k) — > (I,K) [71] given by (A. 6). 

B.1.3 Initial-Initial emitters 

The initial-initial gluon-gluon-gluon antenna is obtained by crossing symmetry from the 
corresponding initial-final antenna function (B.4), with the replacements S12 — > {pi + P2) 2 , 
S13 -> (Pi ~ P3) 2 , s 2 3 ->• (P2 ~ P3) 2 and Q 2 = si 2 + s i3 + s 2 3- It reads [71], 

F°M 3 2 1- 1 ( 8S ^ 4- 8S? 
^3 Kh^g^g) ~ ^TnTvl ~ h 



3 | ° 6 23 


+ 8s 2 3 


+ 


8* 2 2 4 


8s 2 


! S13 






•513 


«23 


^3 + 


4 «23 


+ 


4s?2 




S12S23 


S12S13 




Sl3«23 





2(g 2 ) 2 ^23 812 

+ 12fl3£23 12g 12 g 2 3 12^13^12 
S12 S13 S23 

+24s 12 + 24s 13 + 24s 2 3 + 0(e) , (B.7) 

where the hat identifies the gluons crossed to the initial state. Only the final state gluon 
j may be soft, and it can also be collinear with the initial state gluons i or k. Having 
well defined hard radiators, F®(i,j,k) does not need to be further decomposed. The full 
antenna can be used with a single initial-initial mapping, k) — > (I,K) [71] of the type 
given in Eq. (A. 8). 

B.2 Integrated tree-level three-parton antennae 

In this subsection we give the expressions for the integrated forms of the antennae in (B.l), 
(B.4) and (B.7). The integrated three-parton antennae contain explicit e-poles from the 
integration over the antenna phase space for one unresolved emission and finite remainders. 
They appear in the real virtual channel as single analytic integration of subtractions in the 
(m + 2)-parton channel, corresponding to the double real emission [75], as well as in the 
form of genuine subtraction terms to compensate for oversubtracted poles as discussed in 
Section 2.2.2. Note that the full e dependence in is retained during integration over the 
antenna phase space. 

B.2.1 Final-Final emitters 

For final-final kinematics the integrated antenna was computed in [21] and reads, 

73 

•F 3 °(*i23) = -6lg (6, Sl23 ) + !± + O(e) , (B.8) 
where the colour-ordered infrared singularity operator I^j was defined in Eq. (4.13). 
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B.2.2 Initial-Final emitters 

The full set of integrated initial-final three-parton tree-level antennae were computed in [71]. 
The pure gluon antenna reads, 

^i, j k(sik,xi,x 2 ) = 6{l-x 2 )^6(l-x 1 ) ^-4I«(e,s Ife ) + ^- ivr 2 ^) 



+ 



'Ik 



7^(^)-y^o(xi) + 2P 1 (x 1 )-^- 



- 2(l-2x 1+ xf-x3) _ 2(l-xi + x 2 ) 2 x + ^ 

Xl Xi(l — Xl) 

(B.9) 

The splitting kernels, p® (x), and distributions T> n (x\) that appear above are defined in 
Eqs. (D.l), and (D.4) respectively. The function H(mi, . . . ,m w ;y) denotes the harmonic 
poly logarithms and their notation is also described in Section D. 

B.2.3 Initial-Initial emitters 

The full set of integrated initial-initial three-parton tree-level antennae were computed 
in [71]. For this kinematic configuration the J 7 ® antenna reads, 



^12j(Sl2) = W " 31 Wl " «12) + \ 

+ 5(1 -xi) 



£12 
p? 



\vgg{x2) 



1 - xi 



_ (4 ~A + 2^2 - 1) (log(2) - H(l, x 2 ) - H(-l; x 2 )) \ 

X2 ) 
+ (-1 + X 2 - 2 + 1) P (X1) + ^0(X1)PO(3 2 ) - 2(1 + ,j (1 + , 2) 

| 2(xj-xf + 2)(l + xix 2 ) 2 | 4-xi + 2 + j&TTfc--k 

(xi+x 2 ) 4 ' 1 - x 2 

_ 2 (xf + x\ + (xf + 3) xi + 2) (1 + x x x 2 ) 2 
(l-x 2 )(xi + x 2 ) 4 
3 (xf + xi) (x\ + x 2 ) (1 + xix 2 ) 
(xi +x 2 ) 4 

2 (x 2 x 4 + (xi + l) 2 x 2 + x 2 ,) (1 + xix 2 ) 2x1 (xf + 3) (1 + x x x 2 ) 2 



+ 



+ (xx+x-iY + (x 2 + l)(xi +x 2 ) 4 

+ -. (1 + Xl ^ 2) + ( Xl o x 2 ) + 0(e). (B.10) 

xi(xi + l)x 2 (x 2 + l) 

B.3 One-loop three-parton antennae 

The one-loop antenna functions are obtained from the colour-ordered renormalised one- 
loop three-parton matrix elements according to Eq. (2.59) [21]. These contain explicit 
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poles from the loop integration. The integrated forms for the final-final, initial-final and 
initial-initial cases are available in [21], [72] and [74] respectively. 

B.3.1 Final-Final emitters 

The gluon-gluon antenna functions are derived from the physical matrix elements for H — > 
(partons) [57]. At one-loop, the leading colour term reads [21] 

Fi{l g ,2 g ,3 g ) = 2(lW(e,su) + I^^ 

- f fl(yi2, 2/13) + R(yi3, 2/23) + R(yi2,V23) + y H(0; |yi2|) + y H(0; \yi ?> \) 



11 
"6" 



H(0 ;te l))F 3 »a,2.3) + J- + 5 i- + J- + 5 i- (B.U, 



with yij = Sij/s\23 and where the colour-ordered infrared singularity operator I g ]j was 
defined in (4.13). The function R(y,z) is defined by 

J2 



7T 



R{y,z) = H(0; |y|)H(0; \z\) + — - S(y < O)0(z < 0)vr 2 



+6(0<j/< l)H(l,0;j/) 



+0(y < 0) 



+9(y > 1) 



H 1,0; 



1-2/ 



H 1,1; 



H 1,0; 



y 

+9(0 < z < l)H(l,0;z) 
+0(z < 0) 

+e(z > 1) 



H 0,0; 



7T 
"3 



H 1,0;- 



1-z 



-H 1,1;- 



z 



H 1,0; 



1 



H 0,0; 



1 



7T 
"3 



1-2 

2 



(B.12) 



For the final-final case Sij = (pi +Pj) 2 and the antenna is evaluated in the region where all 
Sij are positive. 

The numerical implementation requires the partonic emissions to be ordered mean- 
ing that the two hard radiator partons defining the antenna are uniquely identified. For 
Fg(l,2,3) the decomposition into subantennae follows the same pattern of the tree-level 
antenna since it can be written as a function proportional to its tree-level counterpart plus 
a function which is not singular in any unresolved limit 

Fi(l, 2, 3) = ^(1, 3, 2) + fi(3, 2, 1) + ^(2, 1, 3). (B.13) 

The subantenna reads 

$(l g , 2 g , 3 g ) = 2(1$ (e, s 12 ) + I«(e, a 13 ) + I$(e, s 23 ) - 2I«(e, Sl23 )) / 3 °(1, 2, 3) 

R(yi2, yw) + R(yi3, 2/23) + R(yi2, 2/23) + y H (°; I2/12 1) + y H(0; |yi 3 |) 



+ H H (0; \ m \)) 2, 3) + -L + J- + 

6 J 6S12 6s 23 9si23 



(B.14) 
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B.3.2 Initial-Final emitters 

The initial-final one-loop three-gluon antenna function can be obtained from its final-final 
counterpart (B.ll) by the appropriate crossing of one of the particles from the final to the 
initial state, i.e. by making the replacements, S23 — > (P2 + P3) 2 > 0, s\ 2 — > (pi — P2) 2 < 0, 
«13 -> (pi - Ps) 2 < and S123 ->■ q 2 = S12 + S13 + ^23 < 0. 

Some care needs to be taken in the continuation, since the final-final antenna function 
is renormalised at p 2 = q 2 = sm, while the initial-final antenna function is renormalised at 
fj, 2 = —q 2 = — S123 [72]. With this in mind, and defining again the j/jj = Sij/si 23 , Eq. (B.ll) 
can be taken over to the initial-final case as: 

Fi(l g ,2 g ,3 g ) = 2(lW(e,sn)+I$!{e,8 13 ) + I^ 

R{yi2,yu) + R(yi3,y2z) + #(2/12,2/23) + y H (°; 1 2/12 1 ) + y H (°; 1 2/13 1 ) 

+ ^H(0; \y 23 \)) F 3 °(i, 2, 3) + -^ + -^ + -^ + -^. (B.15) 

6 / 3S12 3si 3 3s 2 3 3S123 

Just as in the tree-level case, the one-loop initial-final antenna is decomposed symmet- 
rically into subantennae containing only unresolved limits interpolated by the appropriate 
initial-final phase space mapping, 

Fi(l, 2, 3) = f 3 \l, 2, 3) + /!(1, 3, 2), (B.16) 

fi(i g ,2 g ,3 g ) = 2(l«(e, Sl2 ) + I«^ 

- (R(yi2, 2/13) + R(yi3, 2/23) + R(yi2,y23) + y H(0; \y 12 \) + y H(0; I2/13I) 

+ y H(0; \ m \)) / 3 °(1, 2, 3) + J- + J- + J- . (B.17) 
6 / 3si2 6s 2 3 DS123 

B.3.3 Initial-Initial emitters 

The initial- initial one-loop three-gluon antenna function is also obtained from (B.ll) by 
crossing two of the gluons into the initial state with the replacements S23 — > (P2—P3) 2 < 0, 
S12 -> (Pi + P2) 2 > 0, S13 ->• (pi - P3) 2 < and si 23 ^ Q 2 = s 12 + s 13 + s 23 > 0. The 
expression yields: 

Fi(l g , 2 g , 3 g ) = 2 (l« (e, s 12 ) + 1$ (e, s 13 ) + 1$ (e, s 23 ) - 2I«(e, ai23 ))i#(i,2,S) 

- (r(v 12 , 2/13) + R{yi3, 2/23) + ^(2/12, 2/23) + y H(0; I2/12I) + y H(0; \yi 3 \) 
+ H H (0; \y 23 \)) F 3 °(i\ 2,3) + ^ + -^ + -^ + -^-. (B.18) 

O / OS12 OS13 3S23 3S123 

Since the hard radiators are uniquely identified with the initial-state partons, no further 
decomposition is necessary. 
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C. Large angle soft terms for the double real subtraction term da^ 



NNLO 



As discussed in Sect. 2, it is convenient to rewrite the large angle soft terms in a form 
where the first mapping is of final-final type (A. 4). In this case, one should replace the 
wide angle soft contributions given in eqs. (5.9), (5.13) and (5.17) of Ref. [75] by eqs. (C.l), 
(C.2) and (C.3) respectively. 

C.l IIFFFF topology 

ir fa s N\ 2 C(e) 2 , 2 ^ f 

v J K ' " (ijkl) I 

2 ( - s 2i((ii)j) + S 2 i(ii) - Su{( k i)j) + + S m j)i{(ki)j) ~ %)i(fci)) 



+ 2 

x/3^/) ff A,(*i) g )^4(i s ^ 
+ 2 ( ~ s (ii)i«ki)j) + + Sn m)j) - 5i, (fc ,)) 

x/a^i^i^^K^,^,^,^)^^^,^) 

+2 ( ~ 5 («)j'((*j)o + s (ij)j(kj) + s ij(mi) ~ S ^m) 

+ 2 ( ~ S(ik)k((jk)i) + S(ik)k(jk) + S 2 k((jk)i) - S 2 k(jk^j 

xfl(2 g ,i g , Uk) g )Al(l g , %, (Ukji) g , W) g )4 2) (PWi,W) 

+ 2 ( ~ S (H)i((ji)k) + S (li)i(ji) + S 2i((ji)k) - S 2i(ji)) 

x/ 3 °(2 9 , fc fl , 2 9 , «j$fc) s , {U) g )J { 2 \^p^)) 

+ 2 ( ~ 5 i« + 5li2 ~ S 2i(ji) + - Si^a) + 

xF 3 °(i 9 , 2 9 )^(I 9 , § 9 , (j*) 9 , (iT) fl ) J 2 (2) (p^.P^) 

Xi^^K^ (C.l) 

C.2 IFIFFF topology 

,. eFs sr fa s N\ 2 C(e) 2 1JR . x 2 f 

dt7 7ViVLO = -A/iO ( -^JT ) ^7^2 d$ 4(p3,-..,P6;Pl,P2)^ ^ 1 



- 55 - 



{~ S 2j(kj) + S 2j(kj) ~ S Tj(jf) + S 1jCjO ~ ^J' 1 + S 2jl) ^3 (iff. *fl> 2 g ) 

x^(i fl J fl ,(^) 9 ,^ fl )j( 2 >(p (1 ^,p^) 

x^(i p ,2 fl ,(^ fl ,(y) p )j( 2) (p^,p^) 

+ 2 (^(Tfe) ~ S 2jW) + ^(jfc) ~ S lj(jk) + Stfl - %i) ^3(i<?! Ig, 2g) 

xAl(l g , (7T) fl , 2 fl , (^) g )J 2 (2) (p ( ~ r ^) 
+ 2 {- S 2j((kj)i) + S 2 j(kj) ~ Sij(kj) + Sij(i(kj))) fa ^> 



xA°(I fl , 2 fl , ((*j)O p )^ 2) (P^,P ( -j^j) 



+ 2 (-S2i({ki)j) + S 2 i(ki) ~ Su(ki) + fsilgjg, (kl) g ) 

xAl(i g ,(u) q X(m) q )4 2 \p ( iryP l 



+ o i~ S lj{(kj)l) + SljOj) - ^2j(fcj) + S 2j(l(kj))) /3°(2<y, iff, 



iiiy y {ki)jy 

2 

x^(i s , 2 S , (^) 9 )4 2) b(-j,P(^r ) ) 

+ 2 {- s n{(ki)j) + Su(ki) - S 2 i(ki) + S^j^ki))) /3(2ff,is, (feOff) 

x^(i p> (ir) fl ,2 fl ,(^) ff )jW(p ( ^ )>P( -~ ) )|. (c.2) 

C.3 IFFIFF topology 

,. eZf sr /a s iV\ 2 (7(e) 2 1JK , . 2 f 

dc7 7ViVLO = -A/iO ( "^T J ^2d$ 4 (P3, • • • ,P6-,Pl,P2) ^ < 

- (-S li2 + S m + S 2i (ji) - %^ + -S'liO'i) ~ -^3° (iff, 2ff, iff) 

2 {~Sij2 + S Vj 2 + S 2jiij) - S^jTfi + S lj(ij) - S r .^j F 3 °(i ff , 2 g , kg) 
xA° 4 (l g , W) g , (i0 9 ) 4 2) (P&yPvd 

2 (-5ik2 + % 2 + S 2 k(kQ - + 5i fc(w) - -%(^) F 3 °(i 9 , 2 ff , i 9 ) 

x^(i 9 ,^) 9 j 9 ,(ir) 9 )j( 2 )(p~,p~) 



2 



-5ij2 + 5"i[2 + S 2l(lk) ~ + 5 U (, fc ) - S'ji^j) ^3°(iff, 2 9 , 
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(C.3) 



D. Splitting kernels 

In this Appendix, we give explicit forms for the four-dimensional space-like splitting kernels 
used in the paper. Note that we systematically extract a factor of N from the splitting 
kernels, and furthermore, retain only the leading colour contribution. The relevant gluonic 
splitting kernels for the purpose of this paper read [115], 



p° 99 (y) 



p 99 ®p 99 ) (y) 



(D.l) 



^5(1 - y) + 2V (y) + - - 4 + 2y - 2y 2 , 
6 y 

^§~^f) S(l-y)+ 1 j-(2V (y) + ^-A + 2y-2y 2 

+8V 1 (y) - 4 " (0;j/) + 12 - 12y + ^y 2 - 12y H(0; y) + Ay 2 H(0; y) 
1-y 3 



44 4H(0;y) 
3y y 



(D.2) 



t\M = (| + 3 ^) 8(1 ~ y) + 27 + (1 + y) (H H (0; y) + 8H(0, 0; y) - ^ ) 



,3 
+2 



^ - J - 2 - y - y 2 ^j (H(0, 0; y) - 2H(-1, 0; y) - ( 2 

44 



-f Q-y 2 )-12H(0;y)-yy 2 H(0;y) 
+ (2P (y) + --4 + 2y-2y 2 

V y 
x (li " C2 + H(0, 0; y) + 2H(1 ' 0; y) + 2H(0, 1; y) 



where we have introduced the distributions, 



V n {y) 



ln n (l - y) 



(D.3) 



(D.4) 



and adopted the following notation for the harmonic polylogarithms H(mi, m w ; y), rrij = 
0, ±1; the lowest-weight (iw = 1) functions H(m;y) are given by 



H(0;y) = lny , H(±l; y) = T ln(l T y) ■ 
The higher-weight (w ^ 2) functions are recursively defined as 

if m 1; ...,m w = 0, ... ,0 



(D.5) 



H(mi,...,m w ;y) = < 



-ln w y , 



n 



dz/ mi (z)H(m2,...,m„;z) , otherwise 



with 



My) = 



f±i(y) = 



iTy 



(D.6) 



(D.7) 
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We use the CHAPLIN Fortran library [116] to evaluate the harmonic polylogarithms up to 
weight four numerically. CHAPLIN is based on a reduction of harmonic polylogarithms to a 
minimal set of basis functions that are computed numerically using series expansions and 
provide fast and reliable numerical results. 
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